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THE PROPAGATION OF PLANE ACOUSTIC WAVES IN A RADIATING GAS 

By Barrett Stone Baldwin, Jr. 


SUMMARY 

A study is made oj the interaction oj thermal 
radiation and fluid flow in the acoustic approxima- 
tion. An earlier investigation is reviewed wherein 
the disturbances produced in a semi-infinite expanse 
oj radiating gas by sinusoidal motion or temperature 
variation oj a plane radiating wall are analyzed . 
The one-dimensional unsteady-flow equations appli- 
cable to this problem are generalized to include the 
effects of a frequency-dependent radiation-absorption 
coefficient. It is found that a single integro - 
differential equation oj the same form as that 
previously given is obtained. It is demonstrated 
that by a redefinition of parameters , a previously 
given solution applies for a frequency -dependent 
absorption coefficient as well as for a grey gas. 
The solution appears, in general , as the sum of two 
types of harmonic traveling waves: (1) an essentially 
classical sound wave, but with slightly altered speed 
and a small amount of damping, and (2) a radiation- 
induced wave whose speed and damping may be 
either large or small, depending on the frequency of 
oscillation and the condition of the gas. 

The previously given solution for sinusoidal 
motion of the wall is used in superposition to find 
the response to an impulsively moving wall . The 
general results are given in the form of integral 
expressions for the linearized disturbance quantities. 
A closed-form approximation for the velocity dis- 
turbance is obtained. It is found that the main part 
of the response to an impulsive motion of the wall 
propagates initially at the isentropic speed of sound. 
As it progresses, however, the wave becomes dispersed 
and its main part travels at a somewhat slower speed . 
Eventually , the main part of the disturbance shifts 
back to the isentropic speed. Some components of 
the response, associated with the radiation-induced 
wave system, travel at speeds up to the velocity of 


light. As a result, there is a small precursor to the 
main wave front. This part of the response dies out 
exponentially with distance ahead of the main 
disturbance. 

The validity of several approximations used in the 
derivations is investigated. The second-order equa- 
tions are derived to help establish the conditions under 
which the linearized res ults are a first approximation 
to the original nonlinear equations . An integro- 
exponential f unction appearing in the basic equa- 
tions was approximated by an exponential through- 
out the present work. This expedient was introduced 
in early solutions of astrophysical problems. The 
validity of this approach for the present problem is 
investigated by means of a two-term approximation 
and by considering the properties of a n exact solution 
for a grey gas. It is found that the higher derivatives 
of the flow quantities are not given correctly by the 
approximation in a small region near the wall. 
However, valid results are obtained for the flow 
quantities themselves everywhere. 

The integral expressions representing the response 
to an impulsive motion of the wall are evaluated 
exactly in closed' form for a number of special 
cases. This is done to provide information on the 
pressure and temperature fields and as a check on 
the approximate evaluation of the velocity distur- 
bance. Additional checks are made by means of 
numerical evaluations of the integral expression 
for the velocity. 

Finally , the problem of the response to a sudden 
change in the temperature of a fixed wall is con- 
sidered briefly. In this case , a small velocity dis- 
turbance builds up and then decays to zero. Near 
the wall , the gas temperature relaxes exponentially 
from its original value to that of the wall. At 
large distances from the wall, the temperature varia- 
tion is of a type characterized by a diffusion process. 

1 
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INTRODUCTION 

The present work is an extension of the investi- 
gation reported in reference 1 concerning the 
interaction between thermal radiation and fluid 
flow in the acoustic approximation. Tins field 
of study was initiated in 1851 by Stokes (see ref. 
2). His purpose was to show that thermal radia- 
tion does not affect the propagation of sound under 
ordinary conditions. This accomplished, such 
an interaction was not further considered until 
recently. Since 1956 a number of papers on the 
propagation of weak disturbances in a radiating 
gas have appeared. 

Stokes’ investigation was based on an approxi- 
mation appropriate to highly transparent, low- 
temperaturc air. Heat exchange was assumed to 
take place between each element of gas and a 
reservoir at the temperature of the undisturbed 
gas. (For a brief outline of this work see also 
ref. 3.) A parameter that is a measure of the 
rate of heat transfer at a given temperature 
difference has been called the Stokes coefficient 
by later authors. The acoustic equation result ing 
from the inclusion of such a process is a, third- 
order partial differential equation. In the past 
decade two comprehensive surveys on the pro- 
pagation of sound in fluids have appeared in the 
literature (see refs. 4 and 5). In these the treat- 
ment of the effects of radiative heat transfer is 
based on the Stokes approximation. In reference 
C a method for evaluating the Stokes coefficient 
is developed which yields useful information. 
This work utilizes the correct integral expression 
for the radiative heat transfer between elements 
of the gas. The solution is assumed to he of the 
same form as that resulting from the Stokes 
approximation and this leads to an evaluation 
of the Stokes coefficient. In addition acoustic 
wave speeds and damping constants are found 
which compare favorably with the results from 
more recent investigations. 

So far as the author is aware, the first complete 
theory for the effect of radiative heat transfer 
on the propagation of sound far from any obstacle 
was developed by V. A. Prokofiev (see ref. 7). 
Earlier work in this field by Prokofiev is referred 
to in reference 7. This author considers the 
problem of thermally radiating acoustic waves 
in great generality, including the effects of 
viscosity and thermal conductivity as well as 


the smaller effects (for aerodynamic purposes, 
at least) of radiation scattering, radiation pressure, 
and the direct contribution of radiation to internal 
energy. 

In all of the foregoing works, except reference 1, 
attention is confined to the propagation of sound 
waves in the gas far from any boundaries. In 
reference 1 an infinite, plane radiating wall is 
introduced and taken to be the source of the 
disturbances. Following a pattern established in 
previous investigations of chemical and vibrational 
relaxation effects, the influence of the non equilib- 
rium process under study is isolated by neglecting 
other complications. In references 8-11, for 
example, a single chemical or vibrational non- 
equilibrium process is introduced into the govern- 
ing equations with all other processes taken to be 
in equilibrium. Also, in these works attention is 
confined to one-dimensional unsteady or two- 
dimensional steady flows with simple boundary 
conditions. Reference 1 and the present work, 
taken together, represent an attempt to include 
fiow r s with radiative heat transfer in this category 
ol flow fields involving a single nonoquilibrium 
process. In the case of chemical and vibrational 
nonequilibrium, it was found that the same 
governing differential equation applies for either 
process in the small-disturbance approximation. 
It will be seen that the nonoquilibrium effect due 
to radiation does not fall in the same class, 
although there are certain similarities. 

The present investigation follows reference 1 
in assuming that nonequilibriiim effects from 
processes such as molecular transport, dissociation, 
vibration, etc., arc negligible. Radiation scatter- 
ing, radiation pressure, and the contribution of 
radiation to internal energy are also neglected. 
For simplicity, the gas is assumed to be perfect. 
Tiie radiative effects are taken into account on the 
basis of the usual quasi-equilibrium hypothesis, 
wherein a Boltzmann distribution of excited states 
is assumed. The geometrical configuration to be 
considered is that of a semi-infinite expanse of 
radiating gas on one side of an infinite, plane 
radiating wall. Initially the gas is assumed to be 
at rest, in a uniform state, and at a temperature 
equal to that of the wall. One-dimensional 
disturbances can then be produced in the gas by 
moving the wall at constant temperature, or by 
varying the temperature of a fixed wall, or both. 
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In its one-dimensional character, the present 
problem is related to the classical astrophysic-al 
problem of the plane-parallel stellar atmosphere 
(see, for example, ref. 12 or 13). For that case, 
however, the fluid motion is negligible and there 
is no wall. The treatment of radiation in the 
plane-parallel case has recently been extended to 
include the effects of fluid motion and solid 
boundaries (sec refs. 14, 15, and 16). In reference 
1 the wall boundary condition is generalized, and 
the basic equations are used to derive an integro- 
differential acoustic equation for a grey gas 
(absorption coefficient independent of optical 
frequency). This analysis is generalized slightly 
in Section I of the present work by including the 
effect of a frequency-dependent absorption coeffi- 
cient. The resulting linear intcgro-differential 
equation is of the same form as that which applies 
for a grev gas. The only difference is that an 
in tegro- exponential function, appearing as an 
attenuation factor in the equation for a grey gas, 
is replaced by a more complicated function in- 
volving an integral over optical frequency. It is 
found that the radiative properties of the gas 
enter the equation only in their effect on the form 
of the attenuation factor and the value of a mean 
absorption coefficient, both evaluated at the 
undisturbed condition of the gas. To help estab- 
lish that the present linearization is imbedded in 
a successive approximation procedure, the second- 
order equations are derived in appendix B. 

In reference 1, with the aid of a suitable approxi- 
mation to the attenuation factor appearing in the 
radiation terms, the acoustic equation for a grey 
gas is solved for the case ol a black wall under- 
going sinusoidal variations. A similar procedure 
in Section IT of the present work leads to a solu- 
tion of the same form for a nongrey gas. This 
approximation, wherein an integro-exponential 
function is approximated by an exponential, was 
first used in an early solution of the stellar-atmos- 
phere problem (see ref. 12). The validity of the 
approximation for the problem considered here is 
investigated in appendices C and D. Appendix D 
contains discussion of a procedure for finding an 
exact solution for a grey gas. It is found that in 
the exact solution the higher derivatives of the 
flow quantities must he singular at the wall. 
The approximate results do not reproduce this 


effect; however, the flow quantities themselves 
arc adequately approximated everywhere. 

A large fraction of the present effort goes into 
finding the response of the gas to an impulsive 
motion of a wall at fixed temperature. In sec- 
tion III, the solution for this case is found by 
superposition of the sinusoidal solutions using 
Four ior-t run sform theory. The general results 
are given in the form of integral expressions for 
the linearized disturbance quantities. These are 
evaluated exactly in closed form, for certain 
limiting values of the variables, in appendix E. 
A closed-form approximation for the velocity 
disturbance at all values of the variables is de- 
rived in appendix F, and this is checked by nu- 
merical evaluations in appendix G. 

The results from the entire investigation of the 
response to an impulsive motion of a wall at 
fixed temperature are summarized and discussed 
in Section IV. It is found that the main part of 
the resulting wave propagates initially at the 
isen tropic speed of sound. As it progresses, the 
wave becomes dispersed and travels at a slower 
speed, which, for a high gas temperature, becomes 
the isothermal speed proposed by Newton for 
sound waves. These findings are qualitatively 
similar to those for a gas in chemical or vibrational 
nonequilibrium. (For that case, the initial ve- 
locity is the frozen speed of sound, and the final 
propagation velocity is that corresponding to 
equilibrium.) After reaching the slower speed, 
however, the subsequent behavior of the compres- 
sion wave in the radiating gas differs from that 
associated with the other nonequilibrium proc- 
esses. In particular, the main paid of the disturb- 
ance eventually shifts back to the higher isentropic 
speed. In addition, some components of the 
response travel at speeds up t o the velocity of light. 
As a result, there is a small precursor to the main 
wave front at all of its positions. This pait of 
the response dies out exponentially with distance 
ahead of the main disturbance. 

Finally, the problem of the response to a step 
variation in the temperature of a fixed wall is 
formulated and carried to partial completion in 
appendix IT. In tins case, a small velocity 
disturbance builds up and then decays to zero. 
Near the wall, the temperature relaxes expo- 
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nentially from its initial value to a final value 
equal to that of the wall. Far from the wall, the 
temperature variation is of a type characterized 
by a diffusion process. 

In the present work, emphasis is placed on 
fluid-dynamical effects rather than on the physics 
that goes into the determination of the absorption 
coefficient. When evaluation of the absorption 
coefficient is considered, however, common ground 
is established with another class of aerodynamic 
problems that have received increased attention 
recently. In these problems the motivating 
consideration is the evaluation of the heat transfer 
to a space vehicle during reentry into the earth’s 
atmosphere. For this purpose, early theoretical 
and experimental work on the radiative properties 
of air was carried out at the Avco Everett Re- 
search Laboratory (see ref. 17 and refs, thereof). 
A useful semi empirical theory for the emissivity 
of hot air is also given in reference 18. To assess 
the contribution of radiation to reentry heat 
transfer, the effect of the radiation on the fluid 
motion was at first neglected (see ref. 19). In 
reference 20, however, a simplified problem is 
solved wherein the interaction between the radia- 
tion and fluid flow is taken into account. This 
work includes information on other high-tempera- 
ture properties of air, which would be of interest 
in any attempt to compare the results of the 
present work with experiment. 

A related problem, concerning the effect of 
radiation on shock-wave structure, is considered 
in reference 21, It is found that, for weak waves 
in a gas at high temperature, the structure of the 
shock may be determined by radiative beat 
transfer, with n eligible effects from viscosity and 
thermal heat conduction. It is a matter of future 


interest to investigate the relationship of this 
result to the linearized solution for weak waves 
found in the present work. 

Before concluding the discussion of gas proper- 
ties, three other pertinent investigations should 
be mentioned. Reference 22 contains an excel- 
lent review for aerodynamieists of the basic ideas 
involved in the theory of radiative heat transfer, 
including an application of quantum-mechanical 
theory in the evaluation of the radiative properties 
of oxygen. Reference 23 contains a similar 
analysis for a monatomic gas, and the results arc 
applied in a study of shock-tube flow. Finally, 
the results of the present work nmy find applica- 
tion in investigations of very low frequency waves 
in the atmosphere. For information on the 
radiative properties of air under atmospheric con- 
ditions, see reference 24. 

The small-disturbance in viscid-flow theory that 
has evolved over the years takes account of the 
effects of viscosity and thermal beat conduction 
by replacing shock waves and boundary layers 
with discontinuities. The same procedure is 
used in the present work. Thus the gas im- 
mediately adjacent to the wall, when disturbed, 
may arrive at a temperature different from that 
of the wall by virtue of the presence of an opti- 
cally thin thermal boundary layer. For super- 
sonic or unsteady flow, cumulative nonlinear 
effects appear at large distances from the source 
of a disturbance even in the lowest-order small- 
disturbance approximation. Such effects can be 
taken into account by a straining of the coordinate 
system (see refs. 25 and 26). Presumably (lie 
coordinate-stretching process can be carried out 
here, but the matter will not be investigated at 
this time. 



I. ACOUSTIC EQUATIONS FOR A RADIATING GAS 


The linearized equations for one-dimensional 
unsteady flow of a radiating gas are given in refer- 
ence L The coordinate system, shoving the x 
axis to be perpendicular to the bounding infinite 
plane wall, is depicted in figure 1. The semi- 
infinite expanse of radiating gas lies in the direc- 
tion of the positive x axis from the wall. The 
wall is at x=x w (t) i allowance being made for its 
motion. In the present work, the gas will not 
be assumed grey as in reference 1, but the wall 
will be considered black. In that case the acoustic 
equation of reference 1 prior to the inclusion of 
the grey-gas assumption is 


C ) 2 (p 

w 



(y^l Q 

Po 


( 1 ) 


where <p is a potential in terms of which the 
perturbation velocity, pressure, temperature, and 
density are given by the relations 


dx 


(2) 

, cv 

p=-»>Tt 


(3) 

i /av 

dV\ 

(4) 

R\dt 2 7 

bx 2 J 

dp' a v 

Ht~ Po dz 2 


(5) 


The prime is used to denote a perturbation 
quantity. The subscript 0 denotes the condition 
of the undisturbed gas wherein the gas is at rest, 
in a uniform state, and at a temperature equal 
to that of the wall. For example, p 0 is the gas 
density in the undisturbed state, and p is equal 
to p— po- The symbols a 0 , y, and R have the 
conventional meanings so that a 0 represents the 
isentropic speed of sound in the undisturbed gas, 
7 the ratio of specific heats for a perfect gas, and 
R the gas constant. All symbols are defined in 
appendix A. A derivation of equations (1 )— (5) 
is contained in appendix B. 

The quantity Q appearing in equation (1) is 
the net radiant energy absorbed by the gas per 

GGS944— G3 2 


unit volume and time. In general, Q can ^e 
analyzed in terms of its components Q v for par- 
ticular optical frequencies r, where 

(?= J" Q,dv (6) 

An expression for Q> is derived in reference 1 for 
quite general properties of the wall and the gas. 
The subsequent derivation there is specialized to 
the case of a grey gas. In the generalization to 
the case of a nongrey gas, to be treated here, the 
wall will be considered black. The appropriate 
expression for Q u can be obtained by setting the 
quantity equal to one in equations (20) and (21) 
of reference 1. The result is 

Q v — 2Tra lf [B v (T u )E 2 (y v )-{- f B v (T)E 1 (rj y --rj v )dijr 

Jo 

+ P BXhE x {r!-^d^2B v {T)] (7) 

Equations (6) and (7) have not yet been linear- 
ized. They will be rearranged and linearized to 
obtain an expression for Q in terms of for use in 
equation (1). Equation (7) applies in particular 
to the case of a black wall, T w being the wall 
temperature. The quantity a v is the absorption 
coefficient of the gas at the radiation frequency v 



T w =T W (/} ^Temperature of wail 
Figure 1. — Coordinate system. 
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as defined, for example, in reference 22. The 
function B„(T) is the Planck function 

( 8 ) 


and £"1(2), E 2 (z) are the integro-cxponential 
functions 


E n {z)- 


■ rv 

Jo 


2 e 


n= 1,2 


( 0 ) 


For a discussion of the properties of these func- 
tions including a number of integral relations, see 
reference 27. 

The quantity v) p is the optical depth from the 
wall for radiation of frequency v f and is given by 

t) v — f a,(x)dx (10) 

thif (0 


The temperatures T and T that appear in equa- 
tion (7) are the gas temperatures at the positions 
corresponding to 77 „ and 77 respectively. 

A useful alternative form of equation (7) is the 
relation found from it through an integration by 
parts 


<?,=27r«,^ [B y (T a )- 7?,(r) !,„-«]/?* (*,) 


.,07’ 

V>) dr). 

Jo (IT 

df 

~ d f drj dV 


/"I 00 

+1. 


•} 


( 11 ) 


[* B v (T)a„(]v 


Jo 

f* OO 

(12) 

B,(T)dv 


Jo 


= f a{x)dx 

Jx w m 

(13) 


When a v and r\ v are replaced by a. and rj in equa- 
tion (11), and the result is substituted in equation 
(6), the integration over v can be carried out ex- 
plicitly by an interchange in the order of integra- 
tion. This is so because the only remaining v 
dependence is in the quantity B V (E), given in 
equation (8). 

The grey-gas approximation will not be used in 
the following derivation. In the analysis of a real 
gas, for which a v varies with frequency, the 
integration over v cannot be dispensed with so 
readily. However, it will be seen that the averaged 
quantities a and rj are still convenient variables for 
the more general case. 

As a first step toward removing the v dependence 
in equation (11), we note that, for fixed v and t , 
there is a functional relationship between 17 and r\ v 
so that, by a transformation, 77 „ can be replaced 
with rj as the variable of integration. For this 
purpose, the following relations can be established 
by comparison of equations (10) and (13). 

drf (for fixed v and t) (14) 



It follows that 


It should be noted that the partial derivative 
bT/bij, is used to indicate that t is held fixed, A 
further rearrangement of this expression for Q v 
will be made such that upon substitution in equa- 
tion (6), the integration over v can be carried out. 

In the grey-gas approximation ( a v independent 
of v), an average absorption coefficient a and a 
corresponding optical depth 77 are introduced to 
replace and 77*. The method of averaging is 
arbitrary, depending on the weighting function. 
According to reference 16, the Planck mean 
absorption coefficient and optical depth arc de- 
fined by the relations 


dT dr _ 
bX JI ~dv dv 


(10) 


Equation (11) can then be written as 

f” dB.(f) ~ , 07 ' 

jo (IT dr, ‘ 

r°° dB,{T) „ . df \ 

+ ~ J-2\Vy—r).) drj y 

Jv (IT dr, J 


( 17 ) 
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The integral over v could now be completed 
except for the v dependence in the arguments of 
the E 2 function, Tn fact, this can still be done in a 
small-disturbance approximation. This comes 
about because when equation (17) is linearized, 
the v dependent factors no longer depend on the 
temperature variations in the flow field. Then the 
integration over v can be carried out once and for 
all. " 

In the acoustic approximation Q, and hence Q v , 
is needed only to lowest order in the departure 
from the undisturbed gas condition. Prior to a 
disturbance, the gas is at a uniform temperature 
equal to that of the wall. It follows that 
bT 

B V (T W )-~B{T)\^ and are perturbation quan- 
tities. As a result, to lowest order, the other 
factors in each term of equation (17) can be 
evaluated at the undisturbed condition of the gas 
(subscript 0). To accomplish this, the following 
relationships can he established by power series 
expansion 

ll,(T)=B„iT 0 )E B /J— r'+OCT' 2 ) ( 18 ) 

^(T,)=7?.(r 0 )+ f -^p r;+0(7V) (io) 

where 

T' = T-T q , n= T w -T 0 

Since To is constant, differentiation leads to the 
relations 

dB v {T) dB v {T) <1B V (T 0 ) (T ,, 

dT ~ dT c IT o ^ k } 

dT_dT 

dq drj 

Combination of these relations yields 

JWU - JUT) |.o=^ (T’ a - V |,.o)+0(C 2 ) 

dB.{T) dT dJUTo) dT' , n/rr ,„ 
dT dr, ~ dT 0 d v ^ 3 

The symbol 0(T' 2 ) is intended to include second- 

order terms of the form T f -r— and TJ. The last 

07 } 

t wo expressions can be used to write equation (17) 
as 


Q.= 2tt«, 0 E 2 M(T’ u -T' |,=o) 

f- dB,{T 0 ) „ , dT' 

“Jo 'If o Ei(r> * nv) tin dv 

+ (IT ^jE E&-n.) dn+0(T' 2 )] (20) 


The quantity dB v (T 0 )!dT Q is constant, since T 0 
is constant. Therefore dB v {T^)jdT 0 could be 
taken outside the integrals. Instead, in order to 
consolidate the v dependent factors in each term, 
a VQ , which is also constant, will be taken inside 
the integrals in the next rearrangement of equation 
( 20 ). 

The relation between t) v and y is also needed 
only to lowest order in T'. Thus equation (15) 
can be written as 


y v =\ -°-dv+0(T') 

Jo «0 

Since a v jot 0 is constant, it can be taken outside 
the integral and this becomes 

P 4 + o(n=- n+ocn (21) 

Substitution of equation (21) into (20) and 
rearrangement leads to the expression 


„ n (IB,(T 0 ) „ (<**0 \ {T , TV, N 
Q y =2'jra VQ — jjp — h ■> V J U w~ * b=o) 


dB,(T 0 ) 


0 dl o L«o 


(C 2 ) 

*0 J 07 J 

( 22 ) 


The function 2ira v E 2 ( — z) appearing in 

0 di 0 \a Q / 

each term of equation (22) now contains all of 
the v dependence. The subscript 0 indicates that 
Ibis function is to be evaluated at the undisturbed 
condition of the gas. The integration over v can 
therefore be completed for all time, without con- 
sideration of the temperature variations associated 
with particular problems. To this end, the func- 
tion F(ij) is defined by the relation 
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The constant, C, can be chosen so as to obtain a 
formal similarity between equation (27) to be 
derived here and equation ( 33 ) of reference 1 
(result for a grey gas). This is accomplished by 
imposing the requirement that, when a„ Q is in- 
dependent of v (and hence equal to a 0 ), F(rj) 
will be equal to E 2 (r 7 ); that is, 


CE,( V ) = 2r«o —ffi- EMdv 
or 

( 7=2 iraoj^ d,=2ira 0 ~ J" B,{T 0 )<lv 

To evaluate the last expression, use can be made 
of the demonstration in a standard textbook that 
the Plank radiation law (eq. ( 8 )) leads to Stefan's 
law; that is 

P B v (T)dv=-T* (24) 

Jo 7T 


where a is the Stefan-Boltzmann constant (see, 
for example, ref. 28). By differentiation, it 
follows that 

C= 8 aTya 0 (25) 

Use of this in equation (23) leads to 




«’Q dJi,(T 0 ) 

<xq dT 0 



(26) 


Here the parametric dependence of F(rj) on the 
temperature of the undisturbed gas, T U} is indi- 
cated. Henceforth this will not be done. 

When equation ( 22 ) is substituted into equation 
( 6 ), and the order of integration interchanged, the 
result can be written in a form containing the 
combination CF{z) in each term by the use of 
equation (23). Taking C outside the integrals, 
and replacing it by means of equation (25) finally 
leads to the result 


Q=Sa7> 0 [/’«(7’;-7 7 'i, =0 ) 

(27) 

where terms of order T' 2 are neglected. This 
result is formally identical to equation ( 33 ) of 
reference 1 when the latter is specialized to the 
case of a black wall. Tlic only difference is that 
the E 2 (77) function, resulting from the grey-gas 


approximation in reference 1 , is here replaced by 
the function F(r}) defined in equation (26). The 
definition of 77 , equation (13) is the same as that 
used in reference 1 , where it appears as equation 
( 22 ). 

The foregoing procedure can be extended to 
higher order. It is then necessary to define addi- 
tional functions, similar to F(ri) } winch involve 
integration over v of integrands containing higher 
derivatives of B V (TJ) with respect to T bt as well as 
derivatives of (aja) and ( ajp ) with respect to T, 
also evaluated at T= J\. The second-order equa- 
tions are derived in appendix B. 

Equation (27) is now in a form suitable for 
combination with equations ( 1 ) and (4) to ob- 
tain a single expression for the potential <p, 
but in addition a relation between 77 and x is 
needed. This can be found from the definition of 
77 , given in equation (13). To lowest order, this 
becomes 

y=a 0 [x—x w (t)] (28) 


The earlier relations, to be used in the combina- 
tion, are 



dV 


(y-D 0 

a) 

at* 0 

dx 2 ' 


v 

Po 

dT' _ 

X( 

'5V 

al c>y\ 

fA\ 

dt 

11 \ 

dt 2 

7 dx 2 J 



Equation (28) can be used either to replace 77 
in equation (27) or to replace x in equations (1) 
and (4), retaining 77 as an independent variable. 
If the results are extended to higher order, a 
difficulty associated with the transfer of boundary 
conditions is avoided by the use of 77 as independ- 
ent variable. That approach is used in appendix 
B. However, such complications will be ignored 
here, and x used as a variable instead of 77 . Then 
x w {t) is taken to be a small quantity such that a 
displacement through this distance at any point 
in the field leads to higher order terms which are 
neglected in the first approximation. As discussed 
in reference 1 , this is not a uniformly valid ap- 
proximation because of the infinite derivative 
of the E 2 (z) function at a zero value of its argu- 
ment. In the present impulsive piston problem 
there is a further difficulty arising from the fact 
that x w (t) will become large at large time. It 
can be shown that the result obtained by ne- 
glecting z w (Q is valid to lowest order (for values of 
x not too large) if x is measured from the wall 
rather than from a fixed origin. At large x t a 
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straining of coordinates is required to render the 
results valid as in any acoustic theory. These 
points are further discussed in appendix B. 
Substitution of equation (28) into equation (27), 
with xjt) equal to zero, yields 

Q=8<rTfco{^(«oz) [T' a (t)-T'(t,0)] 
-jV[a 0 (r-;)] dx+f'FMx-x)] —■ d*y 

(29) 


Substitution of this into equation (1) and differen- 
tiation with respect to time leads to a relation 

7\2 rp f 

containing Differentiation of equation (4) 

with respect to x relates the latter quantity to an 
expression in terms of <p so that T r can be climated. 
The result is 


j) /dV 
dt\dt 2 ' 



dV\ S(7 l)erTo<*o r — F(a 0 x) 
drv n Po (. 

dT’ w /dV_a|dV\ 1 
dt + W y dxVr-oJ 


Tills equation, derived for a gas with an arbitrary 
dependence of absorption coefficient on frequency, 
is formally identical to equation (38) of reference 1, 
when the latter is specialized to the case of a black 
wall. The only difference is that the E 2 function 
of reference 1 is here replaced by the function F \ 
defined in equation (26). A result not restricted 
to the case of a black wall can be obtained by 
replacing the E 2 function with the F function in 
equation (38) of reference 1. 

The boundary conditions associated with equa- 
tion (30) are 

(i t , 0) =u w (t) = given function of t (31) 

T' w {t) —given function of t (32) 
<p(t f oo) = finite quantity for all t (33) 


II. APPROXIMATE SOLUTION FOR AN OSCILLATING PISTON 


In reference 1 solutions of equation (30) corre- 
sponding to sinusoidal variations in wall velocity 
and temperature were obtained for the case of a 
grey gas {F(z) =E->(z)) by approximating the E 2 
function with an exponential. The same proce- 
dure can be followed in the general case, when F(z) 
is not equal to E 2 (z), by setting 

F(z) *z me~ nz (34) 


When this is substituted into equation (34), and 
the value of n determined by a least squares fit 
of the result to equation (20), it is found that n 
must satisfy the relation 


n ,h - 4 f 

Jo 4 to an Jo 

o 


(IByjdT 0 

-<ru(l+— * 

7 r \ < 


) 


dv = 0 


(38) 


The constants m and n can be chosen, as in refer- 
ence 1, by making the approximation exact in the 
Ross eland limit of strong absorption and by further 
imposing a least squares fit. To obtain the Rosse- 
land approximation correctly, it is necessary to 
match the first moment as follows: 

f~F( S )sdz- f"me-*zi S =% (35) 

It can be shown that this is the correct criterion 
from equation (29). To do this, first integrate 

d 2 T' 

the integral terms by parts so that appears 

under the integrals, A limiting process with <x 0 
going to infinity then shows the result to be pro- 
portional to a double integral of F. By a partial 
integration the latter quantity can be expressed as 
the first moment of F (see ref. 29). 

Substitution of equation (9) into (26) and that 
into the last equation yields 


m 7r 

n 2 ~4^r 0 



dB,(T 0 ) 

dT 0 




(36) 


With an interchange of the order of integration 
the z and y integrations can be carried out to 
obtain 



<lBr(T 0 ) ldT 0 
a ° - aTI 


(37) 


If a V(} is independent of v, it can be seen that this 
reduces to 


or 71 = 1.562 as obtained in reference 1 for the 
case of a grey gas. In that case equation (37) 

reduces to m=~n 2 } also obtained in reference 1. 

o 

In the general case, when the dependence of a„ 0 
on frequency is known, equation (38) can be 
solved numerically for n. The constant m can 
then be evaluated using equation (37). By rede- 
fining the quantities $ and K appearing in equa- 
tions (57 a, b) of reference 1 as 


and 


P 


.Wo 

CO 


(39) 


jr^ \§(y—\)oTl m 
n 


(40) 


one can show that the results of reference 1 in 
terms of £ and K apply for a real gas when m and 
n are determined from equations (37) and (38). 
The validity of the exponential approximation is 
investigated in appendices C and D where it is 
shown that the resulting solutions lead to valid 
approximations of the physical quantities every- 
where. In a small region near the wall, the 
higher derivatives of the physical quantities are 
singular. The exponential approximation docs 
not reproduce this effect. 


10 
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The results from reference 1 for the case of an 
oscillating wall are the following: 

+ rf - cx K £ u)]'“} <41) 


Re [(c 1 <7 I +c 2 <7 2 )e‘“ < ]= 


«0 


where 


+ 0 +%&)«]■ "Y 

a a =yRT 0 


TM 

" T a 


(42) 


(43) 




- (1 TV (1 — <3 2 — iKp)*+4P 2 (l —iKfi/y) 

2(l—iK0jy) 




(44) 


Equation (44) is the solution of a fourth degree 
characteristic equation which results from sub- 
stituting a complex exponential form into the 
in tegro -differential equation. The derivation of 
these relations is contained in reference 1 and 
is further discussed in appendix C herein. 

Equations (41), (42), and (43) apply when the 
wall velocity u w (t) and the wall perturbation 
temperature T’ w (t) are sinusoidal functions of 
time with radian frequency o) and arbitrary phase. 
When this is the case, equations (42) and (43) 
can be solved for the complex amplitudes Ci and 
C 2 . Equation (41) then represents two damped 
sinusoidal traveling waves. The damping con- 
stants and wave speeds are determined hy the 
complex constants Ci and These are given 
in terms of the radian frequency w, the basic 
physical parameters, and the approximation con- 
stants m and n by equations (39), (40), and (44). 

The variations of the damping constants and 
wave speeds through the whole range of pa- 
rameters are discussed in reference 1. As an 
example of these results, the wave speeds i\ and 
v 2 are plotted in figure 2 as functions of the pa- 
rameter w/nat/io f° r 7 = 7/5 and K=4. In the 
figure, the dimensionless quantities i\ /a Q and 
0.01 v 2 /a 0 are plotted to avoid specifying a Q which 
is the isen tropic speed of sound in the undis- 
turbed gas. It is seen that iy is equal to a Q at 
low frequency and also at high frequency. At 
intermediate frequencies the wave speed ap- 
proaches the isothermal signal velocity a 0 j\ y. 


Since the speed of this wave docs not deviate 
greatly from the isentropic speed of sound, it 
has been referred to as a modificd-classical wave. 
The other wave speed, v 2 , varies between zero 
and the velocity of light (taken to be infinite), 
depending on the frequency. The term of equa- 
tion (41) corresponding to v 2 has been referred 
to in reference 1 as a radiation-induced wave. 
The wave speeds are determined from the values 
of the complex constants and c 2 by the relations 

t?i/oo= — [Imaginary Part (ci)] -1 (45) 

V 2 la 0 = — [Imaginary Part (c 2 )V l (40) 

Figure 2 was plotted from data obtained by 
electronic machine evaluation of equations (44)- 
(46). 



Figure 2. — Wave speeds versus frequency of oscillation. 



III. IMPULSIVE PISTON SOLUTION 


The response of a radiating gas to an impulse 
motion of a wall is governed by the integro- 
partial differential equation (30). With the 
exponential approximation (eq. (34)) of the 
attenuation factor F{z) } this is 


d /SV 2 

U \dt 2 a ° 


5V\ 8(7-1)^ 

dx 2 / BpQ 


a 0 m 


{ 


£ ~ ttorO x 



a|dV\ “ 
7 dxVx- o_ 


r 


/^V Qp £)'V \ 
7 das*/ 


+ 


X 


CD 

e 


noro^ - *) 


a Y|?_2l m ,r\ 

bx\dt~ y dxv J 


(47) 


where is a potential from which the gas per- 
turbation velocity, pressure, temperature, and 
density can be found using equations (2)-(5). 
The boundary conditions to be satisfied are 


g;«,0)=„.(0-{" <;<“> (48) 

r - w -{e (!>m «» 

The initial condition 

<p(t } x)=0 (£< 0) (50) 

is also imposed. 

In the results to be given, the perturbation wall 
temperature will be taken zero for all time (0—0) 
and the wall velocity at t ^> 0 will be considered 
constant. However, the required relations for 
more general cases will be developed up to a point. 

Solutions corresponding to general time-de- 
pendent boundary conditions can be obtained by 
superposition of the oscillating piston results set 
down in the previous section. As they stand, 
equations (41), (42) and (43) apply only when 
w»(0 and T' w (t) are sinusoidal functions of time 
with radian frequency o>. For the general time- 
dependent case, the right side of equation (41) and 
the left sides of equations (42) and (43) should be 
32 


integrated over all values of w from zero to <*> as 
follows: 

00 

u(t,x)=— Re I (c 1 Cie ei ‘- z/aa +c 2 C 2 e e,aT,Bo )e ial dc l } 

y J o 

(51) 

Ref (c 1 C 1 +c 2 C 2 )e i ‘ ,t da>=— u w (t) (52) 

Jo Go 

GXs) c, J '“''•} = X <*» 

It is expedient to make use of Fourier- transform 
theory in evaluating f\ and C 2 from equations (52) 
and (53). To this end, the Fourier transform for 
an arbitrary function f(t, x) is defined by the re- 
lation 

7("> 2 -)=X f f(t,x)e~ iw ‘dt (54) 

5T J-co 

where the bar over the function indicates that it is 
a transform. The inverse relation 

/(*iz )=5 f J(a), x)e iwt dco (55) 

— — CO 

follows (see, for example, ref. 30). This definition 
differs from the conventional one in that the ( — i) 
is usually associated with the inversion integral 
rather than the transform as it is here. The rea- 
son for the present choice is to avoid the necessity 
of changing the sign in the appropriate places in 
all of the results of reference 1. When/(f, x) is 
real, it can be shown that the real part of/(a>, x) 
must be an even function of a? and the imaginary 
part odd. When f(t) is also zero for t<C 0, equa- 
tions (54) and (55) can he written alternatively as 

]^,x)=- f j {t,x)e~ iKt dt (50) 

* Jo 

/(£,a-)=ReJ* j (u,x)e iu ‘dw (57) 

By choosing f(t, x) in these forms to be — u u (t) and 
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comparing equation (52) with (57), it is seen that 
equation (52) corresponds to the relation 


C\Ci-\-c<iC2=^- Wjf(w) 
C7 0 


(58) 


Similarly equation (53) becomes 

+R*)G&;)' 


(3), (4) and (62). The results are 
J) '(t,x)= if (O c > 

= -p Q Rc J* (r ie f .“^+r 2 e^“ r / fl »)e < “V/ w 


(03) 


and 


,1 7^) (g9) 

_ 1 0 




| OjJ/ fl{) 


When Unit) and T f w {t ) are given functions of 
t that are zero for t<^0, their transforms can be 
found according to equation (56) and are given by 
the relations 


+ 


= -T 0 Re 


'/.*[( 


(■X) 

f)Vv 


I ajJ’/ffg 


u B («)=- 
n J 

j u w (t)e~ iut dt 

'a 

(GO) 

+(l +0 (04) 

F;.( «)=i j 

f T' w (t)e-^ l (lt 
0 

(61) 

In the derivation of these relations, differentia- 
tions and integrations with respect to t and x are 


These results, substituted into equations (58) and 
(59), lead to two simultaneous equations which 
can be solved for the quantities C\ and C 2 needed 
in the integral expressions for velocity, tempera- 
ture, and pressure, to be given in equations 
(62) (64). The quantities c u c 2 , and /3 appearing 
in equations (5S) and (59) are the same functions 
of co as in the oscillating piston solutions and arc 
given by equations (39) and (44). 

The integral expression for the velocity, equa- 
tion (51), can also be written as 

v(t } aO=!y- f (tq C l e riwTfn "+c 2 C 2 e C2 " x/n ") “tfJu 

=— Ref (cifV Cl " T/ft °+C2C 2e C2UJr/aQ )e iwl (lo3 (62) 

y Jo 

The second equality follows from t lie symmetry 
of the integrand. 

Similar expressions for the pressure and 
temperature can be found using equations (2), 


taken inside the integral with respect to co. Also 
use is made of the fact that the perturbation 
quantities are zero at /<T. These steps can be 
justified only under certain conditions which 
require discussion. It is necessary that the 
integrals with respect to co be integrable at least 
in the sense of the Cauchy principal value. For 
example, if the functions u w (t) and T' w (t) are such 
that the transforms te^(co) or T' w ( to) are singular at 
any point on the path of integration (the real 
axis in the complex to plane), then these singulari- 
ties must be circumnavigated along infinitesimal 
semicircles either above or below' the real axis in 
the complex to plane. In the present problem, it 
can be shown that the path of integration in 
equations (62), (63) and (64) should pass below 
any singularities on the real axis to insure that the 
perturbation quantities a\ ill be zero for t<^0. The 
demonstration is somewhat complicated by the 
existence of branch points. Equation (44) can he 
written as 




- (co 2 — TTalfoo — iKna^o)) T V (co 2 — — iKna^a^y + 4// 2 a 2 a 2 co(co— iKn<x {i aJy) 


2 co ( co — iKna 0 a 0 !y) 


? 


(65) 


GGSD44 G3 3 
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The quantities Cj and C 2 appearing in the inte- 
grands can be found by solution of equations (5S) 
and (59). Upon replacing the quantity 0 with 
the aid of equation (39), the results are 


C x = 


n<X t U ■)”“«(") +iCa 

V 7 / d~ OJf‘2 / ffio 

( 1 +-) ( ° ) * 2 -(l +-) ( 

\ yJxnai/iQ+wci/ \ y/\i 

/ rjy na^\y 

V 7 ' \na Q Oo I 0)Cj/ fif 0 


T' w (u) 


na^ao \ 

/I CKyff y "4“ <jiC 2 J 

( 66 ) 




Yihhf— " a °f° - )<■,-( 1 +-)(- “? -)<•, 

\ 7 / Vn£* 0 (7 o -j~ wC i / \ 7 / y*a 0 tf u + 02C 2 / 

(W) 


It can be seen from equation (65) that there arc 
six branch points in the expression for c x and c 2 . 
Three of these are in the upper half of the complex 
« plane and one at the origin (see sketch). 



The remaining two are in the lower luilf-planc. 
The quantities c t and c 2 will be single-valued 
functions of o> if a branch cut is introduced con- 
necting the four branch points on or above the 
real axis, and a separate branch cut is introduced 
connecting the two brand) points in the lower 
half plane. In the following, the properties of 
Ci and c 2 in the lower half-plane only are of interest. 
It can be seen that Ci on one side of the branch cut 
in the lower half-plane becomes c 2 on the other 
side and vice versa. From equations (66) and 
(67) it follows that the same is true of and C 2 . 
Since the integrands of equations (62), (63), and 
(64) are unchanged by an interchange of the sub- 
scripts 1 and 2, the branch cut in the lower half 
plane can be dispensed with. In general, there 
may be additional regions of nonanalyticity arising 
from u w (o)) and T' w (o)). For the specific problem 
to be considered, there are not. At this point we 


will specialize to those cases where u w (a j) and 
T' w {u) are analytic in the lower half-plane. Then, 
altogether, one can conclude that the integrands 
of equations (62)-(64) are analytic in the entire 
lower half-plane, but contain singularities and 
branch points on and above the real axis. 

The Cauchy integral theorem can now' be used 
to show' that the path of integration in equations 
(62)— (64) should he taken below any singularities 
on the real axis to insure that the perturbation 
quantities will be zero for /<T. Let us assume 
that the properties of F.^(w) and 7*(w) are such 
that the parts of the integrands excluding e iu}t in 
equations (62)-(64) are zero at infinity. Then for 
i less than zero, each line integral involved can be 
closed with a semicircle at infinity in the lower 
half-plane without adding anything to tlie value 
of tlie integral. The resulting closed contour will 
not enclose any poles so long as the path of 
integration along the real axis is taken to be below 
any singularities on the real axis. In that case 
the value of each of the original line integrals is 
zero for 0. 

We will now specialize to the case of impulsive 
motion of the wall with w'all temperature held 
constant. The boundary conditions for this case 
can be written as 


u w (t)--= 


/'! 

\liin L e 


(t<0) 
(t> 0) 


( 68 ) 


T'„(i) = 0 (all 0 (09) 


Substitution of these relations into equation (50) 
yields 

««(»)=- U lim f (70) 

7T ,/() 

and _ 

T{ «)-0 (71) 


Evaluation of equation (70) leads to the result 


u xv {<J)=- Ulim 

7T e->0 



. u 

l — 

7TW 


(72) 


At this point, it is helpful to note that there is a 
basic similarity in the problem. By the intro- 
duction of a set of dimensionless variables, it is 
possible to present the results in a form which is 
independent of the values of the radiation absorp- 
tion coefficient, a 0) the approximation parameters 
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m, ft, and the gas properties a 0i Jl. An appropriate Using these equations, the expressions for the 
set of similarity variables is the following: velocity, pressure, and temperature become 


T ~y 


U 27 




(78) 


V 7+1 

/ l~2~ 

i'=w/y^pY««oo° 


i = (70) 

Aj=ixv^JL 0=1,2) (77) 


i r ( 79 ) 

Poffof 2 tT,J_„ \C, / V 

+( 1 +f)t' vt ]' , "'T < 80 > 

The wave speed parameters Ci, c 2 , and the wave 
amplitude parameters A : , A 2 are given by the 
relations 



As discussed earlier, the line integrals in equations 
(78) -(SO) should be distorted from the real, axis 
near v=0 so as to pass below the singularity at 
the origin. 

The task before us now is to evaluate the inte- 
grals in equations (78)-(S0). This can be done in 
closed form only for certain limiting cases or by 
approximation. The methods used will be demon- 
strated for two cases in this section, and the re- 
mainder will be treated in appendices E and F. 
One question which arises is whether discontinu- 
ities occur in the flow field as a result of the 
discontinuous wall velocity. Such a discontinu- 
ity might occur at r— 0. It will first be shown 
that this is not the case. For this purpose we 
shall concentrate on the expression for the veloc- 
ity given in equation (78). 


It was demonstrated earlier by means of con- 
tour integration that equation (78) will lead to 
ujU— 0 for r<( 0. Therefore, the quantity c ,FT cun 
be replaced by e ivT —e~ iv ^ since the added term 
will be zero except possibly at r=0. Then equation 
(78) can be replaced by 


€) 

U 


f (A l e c ^+A 2 e c ^)(e Ur -e- lrM )- 

’0 (r<0) 

' f” {A l e e ^+A i e c ^)^ 1 dv (r> 0) 

7T ,/ —c V 

(84) 


It can be seen that the real part of the integrand 
in this expression is an even function of r, and the 
imaginary part is an odd function so that equa- 
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2 

U IT 


VT dv 

V 


(r >0) (85) 

For slightly positive values of r, this ran be 
written 

tt(°+, 0 ^2 ]]m S A" V L fa 

L 7 T r->0+ Jo v 

+- lim f _ Re(.lV’ { -M 2 ^) 5 1I L vr ( j v ( 8 ( 5 ) 

The integral is broken into two parts here to pro- 
mote simplifications in eacli part. The point of 
division is, to some extent, arbitrary. By ex- 
pansion of sin vr in a power series and use of the 
fact that the quantity Re [A x exp (c { v£)-\-A 2 exp 
(cr£)] is bounded, it can be seen that the first term 
in this expression for '1/(0+, %)/U is, at most, pro- 
portional to yr nmltiplied by a convergent series 
and hence vanishes in the limit as t goes to 
zero. Then equation (S6) becomes 


|im r 

0/ 7 T r->0-f- Jl/yr V 

(87) 

If the quantity Re (A l e e "*’\-A s p € 2 p *) is expanded 
for large v, it can be seen that only the lowest 
order terms can contribute to the integral in the 
limit. 

For later use, the wave speed parameters c { and 
c 2y given in equation (81), can be partially ex- 
panded for both large and small v in the forms 


i = — i I l — i 


(y —l)ku 
\x , 7+1 


v 2 — iykv-\~- 


(7-l)W 


0(i/ 3 ) (small v) 


2 * 7 ,7+ IV 

V — tykv + g ) 


y ~\ I V 7 ) (large v) 


C—J¥f\ H- 

V v{y—V'" 1 


i (y — 1 )kv 


This result is obtained by rearrangement of the 
quantity under the inner radical in the form 

(v-—iykv—-- 2 -) +4 ( 7 ^ ') v{v—ik) 

=(v 2 —iykv+ y ^ ' ) +4 (+ 1 ) (7—1 )ikv 

and by expanding about a zero value of the second 
term on the right. Further expansion of equations 
(88) and (89) for large v leads to the relations 


<^rS+'>Q>) 

k finite 

(90) 

V© :+'>(:>) 

k finite 

(91) 

of these expressions 

into equations 


(82) and (83), and expansion of the result for large 
v yields 

^1+0( V) k finite (92) 


.1,-1 +0 
T=0 ( 


k finite 


Use of equations (90)-(93) in equation (87) leads 
to 

? r R 

U 7-bH 7T J\lyf~T 

[>+"©]V* 

It can be seen that the contribution of the term 
containing 0^-^is zero in the limit as r goes to 
zero and the integral can he written as 

u( 0+, {) r / y — 1 \ . "1 

y— -™pl -I 2 --)K\ 


lim 2 f" '««> H * 

T->0-r 7T J \ty T V 


By the same procedure as before it follows that 
the contribution to the integral from this integrand 
in an interval 0 to 1/yT is zero in the limit. Add- 
ing tli is interval, again using the symmetry, and 
expressing the cosine and sine in exponential form* 
the integral can be written as 


( 7 +) (7-1) w 


/ 0( P 2 ) (small v) - =e )** lim (- ’ ) | 

{of+durgcl) < 89 > L 


[e l v(t +V + givtr-V — ^-iVfT-£)__^-ip(r+e | 
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Tf we take the line integral to pass below the 
origin, each exponential term can be evaluated 
separately. This is done by closing the contour 
with a semicircle at infinity in the lower or upper 
half plane, depending on the sign of the argument 
of the exponential, and using the calculus of 
residues. The result is the expression 


w(0+i£) 

(J 



T->0+ 


1 (i±L 

2 \|r+€| 


1=L\ 
w-t\) 


The results arc 




(97) 


e— >0 L Ptw J L \ - / - 


(98) 


To evaluate this, the value of £ must be specified 
relative to t before taking the limit. When this 
is done, it can be seen that 

u (0 + , £) __ / 9 0<r<f 

U~ \l o <f<r { } 


Since tlie expansion for large v used in the deriva- 
tion is valid only for finite k, this result may not 
apply in the limit as k —>° although it docs 
apply at all other values of k. In fact, the same 
result can be obtained for infinite k, but the proof 
requires a more detailed procedure which will not 
be given here. 

We have found that the gas velocity is zero 
at slightly positive values of r except at the wall 
(£^0), where the gas velocity is equal to the wall 
velocity. Corresponding results can be found for 
the perturbation pressure and temperature. 
These results are 


//(0+ ,£) f 0 0<r<£ 

PiflJJ \ 1 0 <f<r 


(95) 


RT'(0 + ,& _ /() 0 <r<£ 

aJJ \ ( 7 — 1)/t 9<^<r 


(96) 


By a similar procedure, it can be shown that 
jumps in velocity, pressure, and temperature 
occur at r = £ for all finite £. For the velocity, the 
quantity to be evaluated is 


lim 

e — >0 


L u u J 

“him (=-*) P 

*->0 \ “ 7 T / , 7-00 


[/l 1 e (c i +0 ^+^l 2 e (c 2 +i) ^] “ 


]i]n r RT'(t+€, Q HT’ 

e— >0 L <*'JJ J 

Equations (97)-(99) indicate that the jumps in 
velocity, pressure, and temperature at r=£ 
decay exponentially with distance from the wall. 
This finding is qualitatively similar to that fora 
gas in chemical or vibrational nonequilibrium in 
the absence of radiation (ref. 31). However, 
in the present ease where the effect is due to radi- 
ative heat transfer, the disturbance is not zero 
ahead of the jump (r<(£) as it is for chemical 
nonequilibrium. This will be demonstrated later. 

Tn reference 10, an expression is derived for the 
velocity field far from the wall in the ease of 
chemical nonequilibrium. 1 The same can be done 
in the present problem by expanding the inte- 
grands in equations (78) (SO) for small values of v 
and demonstrating that the contribution to the 
integral from other values of v is zero. The 
evaluation for this and for other limiting cases is 
given in appendix E. 

To obtain a qualitative view of the over-all 
flow field, cither approximation or machine 
computation must be resorted to. Both methods 
have been used to evaluate the velocity disturb- 
ance for the present problem. An approximate 
closed-form solution is derived in appendix F. 


1 In references 10 and 31, the problem treated Is that of the steady supersonic 
flow past a wedge. For chemical or vibrational nonequilibrium there is a 
direct analogy with the corresponding one-dimensional unsteady flow prob- 
lem. Strictly speaking, such an analogy does not exist in the ease of radiative 
heat transfer because of tbc directional propert ies of the radiation intensity. 



18 TECHNICAL REPORT R-13S — NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


The result is 


^( r ,&=0 ( t < 0 ) 



+crfR(l+t+A 

L 2\X 

+| e ~*[ 1 +pi§] (^>0) (100) 

where 

X=(-y/y—l)kt (101) 

b=kj 2+vW2?+l (102) 

The over-all velocity field will be described in the 
next section. 




IV. RESULTS AND DISCUSSION 


Equation (100) is a closed-form uniformly 
valid approximation for the linearized velocity 
response to an impulsive motion of the 'wall. A 
comparison of this result with numerical evalua- 
tions is made in appendix G. Equation (100) 
provides a good qualitative summary of the 
results of the numerical investigation and will now 
be used for that purpose. An evaluation of this 
equation for an intermediate value of the radiation 
parameter (*=-3.0) is shown in figure 3. The 
value *"3.0 corresponds roughly to the value 
K=4.0 used in figure 2. The ratio of specific 
heats y is taken equal to 7/5, but the results 
would be qualitatively similar for any y between 
1 and 5/3. 

Figure 3 is a plot of gas velocity divided by 
piston velocity as a function of time r and distance 
from the wall £. These results apply at all values 
of the mean radiation absorption coefficient, a 0j 
by virtue of the basic similarity and use of the 
radiation mean free path (aj* 1 ) us the unit of length. 
The dimensionless time, r, and distance, are 


defined in equations (73) and (74)_ which are 
T=A /2/(7-i-l )na fi a i) t and £=\ 2 /('y+ Since 

\2/(y+ l) n i s order one, the radiation mean free 
path is approximately equal to oue in units of £. 
The disturbance velocity is plotted as a function 
of time at a series of fixed positions. The r and £ 
scales are broken at several points to produce a 
better visualization of the entire flow' field. Be- 
tween breaks the scales are linear and would have 
their origins at the intersection of the r and £ axes 
if continued hack to the origin. Notice that the 
scale size is quite different in the separate regions. 

At a point located a small distance from the 
wall, the velocity takes a sudden jump at a time 
equal to that required for a signal traveling at 
the isen tropic sound speed to roach the point. 
In other words, initially, the disturbance is a unit 
step propagating at the isen tropic speed. The 
step dies out exponentially with increasing* dis- 
tance from the wall and is replaced by a smooth 
transition from zero at r— 0 to a value of 1.0 at 
T ~+cc. In this process the center of the dis- 
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Fiouiie 3. — Disturbance velocity response to impulsive motion of piston (« 3.0). 
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turbancc shifts toward the path of an isothermal 
signal (t-'v'tS). Eventually, at a large distance 
from the wall, the center of the disturbance 
shifts hack to the path of an iseti tropic signal at 
r=£. The center of the disturbance is taken to 
be the point where the velocity lias reached % of 
its final value when plotted as a function of r for 
fixed £. Tliis point is indicated in each subplot 
of figure 3 by a heavy vertical line under the 
curve. 

At small and intermediate distances from the 
wall, the response is similar to that for a gas in 
chemical or vibrational noncquilibrium (see refs. 
31 and 35). For that case the disturbance is 
initially a unit step propagating at the frozen 
speed of sound. The step dies out exponentially, 
the wave front becomes dispersed, and its renter 
shifts to the slower equilibrium speed. The sub- 
sequent behavior differs from that for a radiating 
gas in that the disturbance continues to travel at 
the slower speed rather than switching back to 
the starting value. There is a further difference 
between the two cases. For a radiating gas, 
there are no characteristics corresponding to a 
finite velocity. Consequently, a precursor extends 
ahead of the path of an iscn tropic signal at t — f 
(see fig. 3). Xo such effect occurs for chemical or 
vibrational nonequilibrium because there is a 
characteristic corresponding to the frozen sound 
speed, and hence no disturbance ahead of this 
line. 

A qualitative physical explanation can be given 
for the response of a radiating gas to the impulsive 
motion as follows: Tteferring again to figure 3, at 
small distance from (lie wall the wave front is 
compact such that its width is small compared (o 
the radiation mean free path (equal to 1 in units 
of £). As a result, the radiative heat transfer 
within the wave front is negligible. The disturb- 
ance is then governed by the isen tropic condition 
and travels at isentropie speed. As it progresses, 
the wave front becomes dispersed by the small, 
but nonzero, radiative beat -transfer process. 
When the width of the wave becomes comparable 
to the radiation mean free path, heat transfer can 
readily occur within the front. This tends to 
hold the temperature constant, depending on the 
intensity of radiation. At high initial tempera- 
ture of the gas (large k), the temperature is field 
essentially constant within the wave front, and 
the disturbance travels at the isothermal speed. 


At large distances from the wall, the wave 
front becomes so dispresed that the radiative beat 
transfer is impeded as a result of reabsorption 
relatively near the point of emission. Thus, when 
the width of the front becomes large compared 
to the radiation mean free path, the adiabatic 
condition applies, and the disturbance travels at 
the isentropie speed. 

There is another useful qualitative explanation 
for the beha vior depicted in figure 3. These results 
can be related to the properties of the solution for 
sinusoidal motion of the wall, given in figure 2, 
as follows. At small £, the wave front is compact. 
Consequently, at a fixed point near the wall, the 
gas velocity varies rapidly with time. A Fourier 
analysis of such a rapid variation would indicate 
a preponderance of high frequencies. In figure 2 
it can be seen that at high frequencies the speed 
of the modified classical wave i\ is equal to the 
isentropie speed of sound a a . This accounts for 
the initial propagation of the compression wave 
at the isentropie speed. At intermediate values 
of £, where we see in figure 3 that the disturbance 
is partially dispersed, the Fourier analysis would 
show the peak amplitude to be at intermediate 
frequency. Figure 2 indicates an approach to 
the isothermal speed at intermediate frequency, in 
agreement with the shift to a slower speed in 
figure 3. When the wave front is further dis- 
persed, at large £, the resulting low frequencies 
lead to a prediction of the observed return to the 
isentropie speed. Obviously, there are gaps in 
the foregoing explanation, if It is not supported 
by other information. But such a point of view 
may be useful in other problems involving a 
Fourier transform that is difficult to invert. 

An understanding of the response to an impul- 
sive wall motion is not complete without considera- 
tion of the effect of varying the radiation param- 
eter k. This quantity is a measure of t lie intensity 
of radiation. It is defined in equation (70), which 

is k= lfi-y/ — g- ^ (y— 1) trTijyTip^ At small 

values of k } the response should become that 
associated with classical acoustic theory, namely 
a unit step in velocity propagating unchanged at 
the iseti tropic speed. At large values of k, one 
would again expect the classical result, but with 
the isentropie speed replaced with the isothermal 
speed, since the large radiative heat transfer 
associated with small temperature differences 



PROPAGATION' OF PLANE ACOUSTIC WAVES IN A RADIATING GAS 


21 



Figure 4. — Disturbance velocity response to impulsive motion of piston (* = 0.001). 


would liold the temperature constant. These 
expectations are realized in a sense, but there are 
singular perturbation effects at both small and 
large k resulting in a somewhat more complicated 
picture. Such effects could be predicted, at least 
in part, from the results given in reference 1 
(response to sinusoidal motion of the wall). 
There it is shown that, for large k, the classical 
result of undamped sine waves traveling at the 
isothermal speed is obtained, except for very high 
or very low frequencies (see fig. 3 of ref. 1). 

At small Ic, equation (102) indicates that the 
quantity b becomes 1.0, and using equation (101), 
equation (100) can be written as 

“ vt ^ r T — I =!, (1— exp [ — (V"V — !)*€]) 

( erf {^ 4 ^r} +erl {2KV7-1W 2 } ) 

+|exp[-(VT-l)**](l+j^E|) (103) 

At all values of £ except large values comparable 
to the quantity exp l— (yy— 1)££] is equal to 
1.0 so that equation (103) represents a unit step 


propagating at the isen tropic speed (the path of 
an isentropic signal is at 7— £). This is the classi- 
cal result which is expected for k — 0 when radiation 
effects are absent. However, at large values of £ 
comparable to 1 jk t the step dies out exponentially 
and is replaced by the error function variation. 
The distance from the wall at which this occurs 
goes to infinity as k goes to zero. At such dis- 
tances, the transition in velocity still occurs in a 
region narrow compared to the distance, and the 
transition is centered on the path of an isentropic 
signal. This result is plotted in figure 4. Since 
the unit step propagates unchanged until the 
distance from the wall becomes large of order 
k~\ the region where it does change is brought 
into focus by plotting the velocity as a function of 
kr and ££ rather than r and £. Two other features 
in figure 4 require explanation. The width of 
the region where the velocity transition occurs 
depends on the value of k and is narrow compared 
to the distance from the wall. Therefore, a 
specific small value of k equal to 0.001 is used. 
Also the scales are broken in several places with 
different scale sizes in each region. 

For small k we see in figure 4 that the disturb- 
ance remains centered about the path of an isen- 
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Figure 5. — Disturbance velocity response to impulsive motion of piston (k — 1000), 


tropic signal (kr=k^), The radiation merely 
causes a dispersion so that the unit step is eventu- 
ally replaced by a smooth transition. 

At large values of k (high temperature), equa- 
tion (102) indicates that the quantity b is equal 
to k and equation (100) becomes 


figure 3. A value of k — 1 000 was used to compute 
the points, but the graphs would remain similar 
for larger values. There is a boundary layer in 

which £ is small of order i where the appropriate 

inner variables are kr and &£, as shown in figure 5. 



Three ranges of distance from the wall must be Within the boundary layer this result includes a 

considered in this case. All three ranges are step which propagates at isentropic speed and dies 

shown in figure 5. Here the scales are broken in out exponentially. The error function terms 

several places and are linear between breaks as in represent a smooth transition centered about the 
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path of an isothermal signal (fcr=V 7 &i;). When 
becomes large compared to 1.0, the latter varia- 
tion replaces the step completely and we have 


Ua 

U 


k-¥*> 

A£ fixed 
lir fixed 



lcT—^yki \ \ 


(10G) 


Tins function represents a transition from zero to 
one in a narrow region centered about kT=\ ; yk%. 

In the next range of distance from the wall, the 
appropriate variables are r and £. These are taken 

to be large compared to ^ but small compared to 

k. Tn that case exp [— (V7 — !)£/£] = 1 - 0 , CX P 

[-(V 7 -l)i£]= 0 > and erf 


Equation (105) then becomes 


Ua 

V 


2 

£ fixed 
t fixed 



y\ 

2[(V7-1)?] 1/2 // 


(107) 


(104) simplifies to the expression 

=K—p I]} 


CO 

|£/fe fixed 
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1+crf - 

i r t"H l 
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-/ 


r— j- I 

r 

^1-Trr 

"V 


1+erf h 

1 V»- «•/ 1 

1 . £ 5 j 

► 

L Aj 


l 2 [ (V1 " 1> 1] . 

> 


(108) 


At e/i«1.0, this reduces to equation (107), and 
hence represents a unit step propagating at the 
isothermal speed. As £/& increases, the step dies 
out exponentially and is replaced by the error 
function variation indicated in the first term. 
In this process the precursor, extending ahead of 
the line r=£, reappears (see fig. 5). At values 
of £/£>> 1.0, equation (10S) becomes 


By cancellation of yX from the numerator and 
denominator of equation (10G), it is seen that the 
result at the outer edge of the boundary layer 
matches equation (107), which applies outside 
the boundary layer. In terms of the physical 
variables r and £, equation (107) represents a 
unit step propagating at the isothermal speed. 
At least it is a step in the limit as k goes^to^ in- 
finity , since the error function structure then 
becomes compressed into an infinitesimal width. 

Thus, except for very small £ of order ^ and very 


large £ of order k (sec below), the expected classical 
result of a unit step with isothermal speed is 
obtained. Then the step which appears at 
in the middle and outer regions in figure 5 
has an error function structure, not resolved in 
the graphs. 

The last range of distance from the wall, of 
interest for the case of large k, comes into focus 

T t 

in terms of the variables -r and jr When these 


arc taken to he large compared to equation 


r 


ujjill 


U 


£/fc fixi'd 
rjk fixed 
£//;-» < x > 


2 


1+orf 


V. 



J J 


(109) 


This result yields the asymptotic behavior at 
large £/£. The transition of the velocity dis- 
turbance from zero to one starts ahead of the path 
of an isen tropic signal (t/£=£/£) and is com- 
pleted in a width of order [(\7 — l)£/£p ~ t(v7 
— 1) r/£p. Thus the width of the wave front 
grows parabolically at large distances from the 
wall; whereas the separation between the paths of 
isothermal and isen tropic signals increases at a 
faster linear rate. 

Tn retrospect it can he seen that the behavior 
of the disturbance for large k is qualitatively 
similar to that for the intermediate value of 
£ = 3.0, described earlier and pictured in figure 3. 
But for large k } the switch from isentropic toward 
isothermal speed of the wave front occurs at 
smaller distances from the wall, is more complete, 
and persists to a larger distance from the wall. 
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Numerical investigations have not been made 
of the pressure and temperature fields. However, 
the behavior of these quantities can be described 
qualitatively using the results from the limiting 
cases discussed earlier and in appendix E. Since 
the results for large k are similar to those for 
intermediate values, only large k will be considered 
hero. 

Referring to figure 5 , at values of ££<<1.0 (i.e., 
near the wall) the velocity undergoes a unit step 
at £r=££ and thereafter remains constant. The 
dimensionless perturbation pressure p'jptfiqU also 
takes a unit step at this point, but does not remain 
constant thereafter. From equations (F.l)- 
(F. 3 ) and (F. 14 ), it can be seen that in a time of 
order £r=1.0, the pressure at the wall drops to a 
value of 1 /\y. It remains at this value until a 
large time of order r=k, when it slowly returns 
to a value of 1.0 and thereafter remains constant. 

At a point located outside the boundary layer 
such that £~ I «£<<£ the velocity undergoes an 
error function variation from zero to one in a 
narrow region near r=V 7?. The perturbation 
pressure follows a similar variation from zero to 
l/v'7 in this region. But rather than remaining 
constant thereafter as does the velocity, the 
perturbation pressure eventually rises to 1.0 at a 
time of order r—k . 

At large distance from the wall, where £ is 
large compared to k, the pressure variation follows 
that of the velocity. This is an error function 
variation from zero to one in a relatively narrow 7, 
region near r=£, and no further change. In all 
three regions the final value of the perturbation 
pressure is 1.0, which is the same as that which 
would occur in the absence of radiation (£— 0). 

Equations (F.1)-(F.3) and (F. 15 ) can be used to 
find the behavior of the temperature for the case of 
large k. At the wall the dimensionless perturba- 
tion temperature RT'ja Q U undergoes a step at 
t = 0 similar to that taken by the velocity and 
pressure, but the amplitude is equal to (7-I )/y. 


This is the same variation as would occur in the 
absense of radiation. ITow r ever, in a time of 
order kr~\. 0, the perturbation temperature drops 
asymptotically to zero by radiation to the wall, 
which is held at constant temperature. It should 
perhaps be reiterated that, in the present inviscid 
approximation, the temperature of the gas adjacent 
to the wall need not be equal to the wall tempera- 
ture at all times because of the presence of an 
optically thin thermal boundary layer. 

At a fixed value of £ of order one, the pertur- 
bation temperature (as a function of r) remains 
zero in the neighborhood of r=< y £ where the 
velocity and pressure rise. It subsequently rises 
slowly and falls again to zero at some time greater 
than r~k. 

At large distances from the wall (£>>£), the 
temperature participates, along with the pressure 
and velocity, in an isen tropic variation near 
r=£. Thus the perturbation temperature, UT'l 
OqU, rises to a value of (7— 1)/7. In the subse- 
quent time, while the pressure and velocity re- 
main constant, the temperature slowly returns 
to its initial value of zero in a time greater than 
r=k. The final perturbation temperature in all 
three regions is zero owing to radiation to the wall, 
which is held at constant temperature. 

In reference 1, a solution corresponding to 
sinusoidal temperature variations of a fixed waJI 
is discussed, as w'ell as the results for sinusoidal 
motion of a wall at constant temperature. The 
problem with impulsive time dependence, corres- 
ponding to the former, would be that of an im- 
pulsive temperature variation of a fixed wall. 
This is an interesting case. Although no complete 
study of the problem has been made, a qualitative 
description of the results to be expected is given in 
appendix II. 

Amiss Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Aug . f ), 1962 



APPENDIX A 

TABLE OF SYMBOLS 


a 


«o 


ttj, fl 2 


A 


A h A a 


b 


bit b-2, b$ 


B 

By 

C 

Cl, C 2 


C 

C, 

erf (jj) 


E n {ri) 


f(>, x) 
/(«;*) 
m 

F(n) 

< 7(0 


quantity equal to (£j3/7)/[l+20 2 / 

(t+1)] 

isen tropic speed of sound in undis- 
turbed gas 

positive finite nonzero constants 
defined in equations (E.S) and 
(E.0) 

complex constant defined in equa- 
tion (C.3) 

complex wave amplitudes defined 
in equations (77) and (78) 
imaginary part of one root of 
v 2 —ik+ 1=0 (sec eq. (F.19)) 
also used as a parameter equal 
to aoCtoi co in appendix C 
positive finite nonzero constants 
defined in equations (E.S), (E.9), 
and (E.13) 

complex constant defined in equa- 
tion (C.4) 

Planck function, see equation (8) 
velocity of light 

complex constants containing 
wave speeds and damping con- 
stants of m odified -classic al wave 
and rad i at ion-induced wave, 
respectively 

constant equal to SaIUar 0 
complex wave amplitudes 
error function equal to (2/yV) 

P e~ s2 ds 
J 0 

integro-cxponcntial function equal 

to f 
Jo 

arbitrary function 
Fourier transform of/(f, x) 
dimensionless wall velocity defined 
in equation (B.S) 
attenuation function defined in 
equation (26) 

d i m cnsionless perturbation wall 
temperature defined in equation 

(B.0) 


h 

Planck’s constant, also enthalpy 
per unit mass 

im 

complex quantity defined in equa- 
tion (C.2) 

Im 

imaginary part 

h, h, h 

integrals defined in equation (E.7) 

k 

Boltzmann’s constant, also used 
as a quantity equal to 16 

V (7 + 1 ) /2 (min) (7—l)<rTg/ 

(7 Rp 0 a Q ) 

K 

equal to 1 6 (m/n) (7 — l)<r 73 / 
( 7 ip a a 0 ) 

m, f n 

quantities defined in equations 
(34), (37), and (38) 

0(e) 

quantity ’of order e 

V 

gas pressure 

p' 

perturbation pressure 

po, Pi, Pi 

components of pressure in expan- 
sion (see appendix B) 

Q 

net heat absorbed per unit volume 
and time due to radiation 

Qy 

net heat absorbed per unit volume, 
time and frequency interval due 
to radiation 

(QM) 1 , (QMh 

components of ( Q/a ) in expansion 
(sec appendix B) 

r 

variable equal to kv 

B 

gas constant 

Re 

real part of 

t 

time 

T 

gas temperature 

To 

undisturbed gas temperature 

T' 

perturbation temperature 

T h T 2 

components of temperature in 

expansion (see appendix B) 

T 

perturbation temperature at x—x 


wall temperature 

m co) 

Fourier transform of perturbation 
wall temperature 

u, u' 

gas velocity 

U U u 2 

components of gas velocity in 
expansion (see appendix B) 

ujt ) 

wall velocity 
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U A 

«.v 

V 

x 

x„ 

X 

a 


« o 


jS 

T 

c 

V 

V 




Fourier transform of wall velocity VvQ) Vn 

approximate gas velocity in im- 
pulsive piston problem # 

remaining part of gas velocity 


after subtraction of approximate © 

part (evaluated numerically) K 

constant wall velocity at £]>0 in 
impulsive-piston problem v 

Cartesian coordinate 
coordinate at the position of the £ 

wall P 

variable equal to (Vy— 1 )&£ Pu 

Planck mean radiation absorption Pn P 2 


coefficient defined in equation 
(12) o- 

value of a in the undisturbed gas 
components of a in expansion (see 
appendix R) 7 

frequency dependent radiation ab- ^ 
sorption coefficient w 

parameter equal to natfL 0 /co 
ratio of specific heats 
a number small compared to one 

C X A 0 

variable equal to adx 

J*w<o 1, 2 

dummy variable of integration 

fn a / 

variable equal to — drj 

Jo a 


components of 7j y in expansion 
(see eqs. (B.14) and (B.15)) 
variable equal to t (used to 
facilitate transformation) 
value of wall temperature at 0 
parameter equal to 8 ( 7 — 1 ) a Tlf 

frequency of electromagnetic radi- 
ation 

variable equal to \' 2 /(y+ l)na {j x 
gas density 

density of undisturbed gas 
components of density in expan- 
sion (see appendix B) 

Stefan Boltzmann constant equal 
to 5.673 ><10“ 5 erg cm -2 deg -4 
sec -1 

variable equal to y' 2 /(y 4 - l)natfi 0 t 
velocity potential 
radian frequency of oscillation in 
oscillating piston problem 

SUBSCRIPTS AND SPECIAL SYMBOLS 

denotes the undisturbed gas con- 
dition (as in p 0 ) 

components of expansion used in 
appendix B, for example, p b p 2 
denotes perturbation quantity (as 
in p' = p— p (J ) 



APPENDIX B 

SECOND-ORDER EQUATIONS 


To help establish tlmt the linearization used in 
the text is imbedded in a rational successive- 
approximation procedure, the second-order equa- 
tions will be derived in this appendix. The small- 
disturbance expansion could be carried out in a 
number of different ways. As stated in the text, 
nonuniformities can be partially avoided by 
transforming to the variables t, r, in place of t, x. 
This procedure will be used here. 

The one-dimensional unsteady inviscid-flow 
equations arc 


dp 

dt 


+ u 


dp dn 
dx +p dl= Q 


(B.l) 


(B.8) 

7 

TM) = T Q +eT 0 <j(t) (B.9) 

where /(0 and <j(t) are given functions of t with 
maximum values of order one. Then e is a 
dimensionless parameter that goes to zero in the 
limit of a vanishingly small disturbance. It 
will be assumed that all quantities can be ex- 
panded in powers of e at least to second order, 
for example, 

T=T Q +eT l + t 2 T 2 + . . . (B.10) 


dn du dp 

p-^r J rP u 3 ^+^=° 


d^ 


dx dr 


dh . dh dp d p 


dt 


dx 


(B.2) 

(B.3) 


(see ref. 1). These are supplemented by equations 
(6)-(17) of the text, which define Q , and the 
equations of state for a perfect gas, which are 


We will first expand the quantity (QJa) to 
second order in e, and return later to the other 
equations. The variable (QJa) is chosen, rather 
than Q , to promote simplifications which will 
become apparent later. Equations (6) and (17) 
can be combined and written as 

5=2, r (f ) { [ T) I, _o]7? 2 (i|,) 


h= y ,RT (B.4) 

7 — J 


J o (IT 






d f 

dr, 


(trj 


The boundary conditions are 

4t,*ic(t)]=ujt) (B.O) 

u(t,x ) finite at x >0 (B.7) 

The wall velocity u w (t) and wall temperature 
T^t) are taken to be given functions of t ( T w (t ) 
appears in the expression for Q). 

We wish to expand the flow quantities about 
a zero value of an appropriate parameter that is 
a measure of the magnitude of the disturbance. 
For this purpose, u w (t) and T w (t) can be expressed 


+ rdjuf) 

Jv df 




df 

dy 


dr^^dv 


(B.l 1) 


where 



(15) 


The integral over v can be carried further by 
expansion of equations (B.l I) and (15) in powers 
of e. For this purpose, the following expansions 
are needed 


°^) +e (^) +0(e 2 ) (B.l 2) 

a \oi/( ) \«/i 

^=^o+ €1 ?n^^( e2 ) 


as 


(B.l 3) 
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By substitution of equation (B.12) into (15), and 
noting that (— ) is constant, it can be established 


’’■■"(I), ” 

--m 


(B.14) 


(B.15) 


Power series expansion of the E« function yields 
E 2 (r) y )=E 2 (vy 0 +eVy l + • • •) 

= ^(^ C )+-^ F °- «7,,+0(e s ) 

The relation dE 2 (v) I dr) — — E x (y) can be used to 
write this expression as 

E-i{i) t ) =E 2 ( v , 0 ) — ^(^jjrj+O^ 2 ) (B.16) 

Since the E x function is logarithmically singular 
at a zero value of its argument, the radius of 
convergence of the power series expansion of the 
E 2 function goes to zero as r\ VQ goes to zero. How- 
ever, this docs not invalidate the expansion, since 
6T}v l goes to zero at the same rate as does the radius 
of convergence. In other words, equation (B.16) 
remains valid at all values of 77. This can be seen 
by noting that as 77 goes to zero, where E x (77^) is 
logarithmically singular, equation (B.15) shows 
that 77^ goes to zero algebraically. Then the 
combination ^1(77^)77^ remains finite and, in fact, 
goes to zero as 77 goes to zero. The same would 
be true of all higher-order terms in equation (B.16). 
Using equations (B.12)-(B.16), the quantity 

— E 2 (t) v ) in the first term of the integrand of 
or 

equation (B.ll) can be written as 

+<{©. a [©/]-(!).■ *■[(!).’] 


we will continue to concentrate on the first term 
of the integrand of equation (B.ll). By power 
series expansion of the function B y (T ), given in 
equation (8), and using equation (B.10), we find 

B r (T) = B,(T 0 ) + e 


2 rdB y (T 0 ) 1 <! 2 B,(T 0 ) 1 

* L dT 0 1 ' ' 2 dT% 'J 


+ ■ ■ ■ (B.18) 


From this and equation (B.9) it follows that 
B.{T W ) - B,{T ) |,_ 0 = e I W) - r,|,-.] 


i f dB y (T 0 ) j, r 1 

\ ITT 2 


1 d 2 B y (T Q ) 

2 dTl 




?U-ol • • • 


(B.19) 


With the assumption of a Boltzmann distribution 
of states, the quantity a v /a is a function only of 
temperature (for a fixed value of v), since it is 
proportional to the mass absorption coefficient. 
This is the reason for choosing to expand (Qfa) 
rather than Q. Expansion of a v fa in powers of 
(T—T q ) and use of equations (B.10) and (B.12) 
leads to the result 


d — (T q ) 

a 


(B.20) 


Then, since T 0 is constant, it follows that 


I'©/;- 


d^(T 0 ) 

a 


(B.21) 


Use of equations (B.12), (B.17), (B.20), and (B.21) 
and a similar expansion of the integral terms of 
equation (B.ll) leads to the results 

Q!a=e(Qj«) l +e 2 (Q/a) 2 + . . . (B.22) 

(<?/«)i=--f ITogiQ-T^Fiv) 

a Q t. 

,i? } <M) 


For the purpose of expressing ( Q/a ) in terms of T, 
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(<?/«) 


J, 


■£S -T 2 \, =0 F(v) 

Oi 0 V. 


-J 

ff’-g m,-5i) 4 

'o \V — V\ 07] J 


+^/fzV(0-z , ;i,- 0 ]F,o») 

« o l 


-J 

r (s=? 

U v—y\ J 




(B.24) 


where the functions F , F r and F n are defined 
by the relations 
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«0 


faA tIB P (Tp) p r /a, 


F <’>“ 2 '.r &X : lfr E 


f / « A 

- \ a / o _ 




dv 

(B.25) 

dv 

(B.26) 


Cn 

«o 


F„(i/)=2x( 


(ITa JT, 




* 


(13.27) 


Tlie constants (7, Cj, and C/ 7 are arbitrary. In 
the text, (7 was chosen such that for a grey gas 
F( rj) becomes equal to E 2 ( 77 ). A similar choice 
here leads to 

<7=Scr7> 0 (25) 

<7,= 12*7% (B.28) 

(7//=8ct7 7 oCKo (B.29) 

Then for a grey gas F i (t]) also becomes equal to 
E 2 ( v ). The function F /7 (i?) introduces the effect of 
temperature dependence of the mass absorption 
coefficient into the equation. In the case of a 
grey gas, this function is zero. 

In the foregoing procedure there are no obvious 
difficulties that would prevent an extension of the 
expansion to higher order in €. It can be shown 
that an attempt to expand E 2 (*?„) using x as the 
basic variable instead of rj will fail because, in 
general, x is not zero when t] v is zero (see eq. (10)). 


A new variable x=x— x w (t) could probably be 
used, but the second-order part of ( QI<x ) would 
contain additional terms depending on the density. 
Altogether, y is probably the most convenient 
variable, when the expansion in e is extended 
beyond the linear approximation. 

It can be seen that the transformed coordinate 
system b 7 ? is 110 1 an inertial reference frame. 
This complicates the expansion of the continuity, 
momentum, and energy equations. However, it 
will be seen that this complication is similar to one 
which occurs in classical acoustic theory also, and 
can be treated. To facilitate the transformation 
procedure, t will be replaced by 6. Then we wish 
to express equations (B.1)-(B.3) in terms of the 
variables 6,7] defined by 

6=t (B.30) 



(13) 


The transformation of derivatives is, in general, 


d __ dfl d . dr; d 
dt~dt dd dt dr; 

5 _b0 5 . dr/ d 
dx dxdfTdx d* 


CB.31) 


(B.32) 


These can be partially evaluated using equations 
(B.30) and (13) to obtain 


d _J).drj _d 
d t dO'dtdr) 


(B.33) 


j> = A 

dx ^ dr; 


(B.34) 


drr 

The quantity ~ can also be evaluated in terms of 
t and x using equation (13). However, in the 
transformed equations, ^ must be expressed as a 


function of 6 and r). For that purpose, equations 
(B.30) and (13) must be inverted to obtain expres- 
sions for t and x in terms of 6 and t\. Since a is a 
function of temperature and density, the inversion 
will depend on the condition of the flow field, in 
general. However, it turns out that with an 
expansion about e=0, the inversion can be found 
to any desired order in e by a “boot strap” type 
procedure. 
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Differentiation of equations (B.30) and (13) 
yields 


dO=dt- 

dt-\-adx 


where 


dt) _ 

dT 


dx w 


r 

J x * 


da(t, x) 


(B.35) 

(B.3G) 

dx (B.37) 


}x w (t) dt 

Inversion of equations (B.35) and (B.3G) leads to 
t^S (B.30) 


dx=~ dr)—- d6 
a a dt 


(B.38) 


Using equations (B.S), (B.30), (13), and (B.3S), 
equation (B.37) can be written 


5 7 ] 
dt 


=- t a w ^M+ P- 

7 ,/o « 


da(t,x) 

dt 


<lh (B.39) 


To express this as a function of 8 and ij, a (runs- 

~&CC 

formation of — according to equation (B.33) is 

needed. But ^ itself appears in that equation. 

The process can be carried out, nevertheless, 

because, for evaluation of ~ to a given order in e . 

dt & 


dt 


is only needed to the next lowest order. For 


this purpose, an expansion of a is required as 
follows 

of=« 0 -(-€Of 1 -f 0(e 2 ) (B.40) 

Substitution in equation (B.33) and using the 
fact that a u is constant leads to 

^ j drj Dai | /t> 4 1 \ 

di =e a*+ € d7dU 0(e) (B41) 

When this is substituted in equation (B.39), it is 

seen that ^7 is of order e. Hence 
ot 


dt =t d0+ O(e } 


(B.42) 


and if follows that 
(Jo 


dr) 

dt~ 


-6«0 ' 


m+ e P 

an Jo 


da x 

~d0 


4+0(e 2 ) (B.43) 


This procedure could be carried to higher order, 
but for the second-order results the derivatives 
need only be transformed to order e , since the 
quantities on which they operate are of order e. 

Expanding the flow quantities in powers of e 
and using equations (B.33), (B.34), and (B.43) 
leads to an expression containing e and e 2 terms. 
Equating the coefficient of the e term to zero and 
replacing Q with t yields 


dpi 

dt 


+«„P„^ = 0 


du l 
Po dt 




d hi dp l 

po At dt a " 



(B.44) 


(B.45) 


(B.40) 


Similarly, equating the coefficient of the e 2 term to 
zero leads to the expressions 
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du, a, dWi 
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P 0 


dll2 dpi 
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(B.47) 


(B.48) 


+: 


«0 


C” da, , A 1/ d/i, dp,' 

J 0 dt <l7> ( P ° d^-fr, , 


dh, a , / 

~ Pi art; v Po 


dhi dp A 
dt dt ) 


(B.49) 


The quantity a t appearing here should he expressed 
in terms of the other variables. This can he done 
by expanding the quantity (a/p) in powers of 
(T— Tq). It can be seen in equation (12) that (a/p) 
is a function of temperature alone, since the mass 
absorption coefficient (a„/p) is a function of tem- 
perature alone in the assumed absence of non- 
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equilibrium processes other than radiation. It is 
found that a x is given by 


d a -(T 0 ) 

Pi i P T 7 

«i=«o rPo — yf 1 1 

p o 


(B.50) 


The expansion of the equations of state is 
straightforward and need not be discussed. 
Finally, using equations (B.6), B.7), and (B.8), 
the boundary conditions can be written as 


/(0 

«2!*-C=0 


(B.51) 

(B.52) 


and all quantities must be finite at > 00 . The 
foregoing expansion procedure could evidently be 
extended to any desired order in e. 

When solutions of the foregoing equations in 
terms of t and ?j are found, a parametric repre- 
sentation would result from also expressing x as a 
function of t and tj. Such a relation follows from an 
integration of equation (B.38) using equation 
(B.43). If 6 is replaced with t and equation (B.8) 
used, the result is 




c(t) — \ f [2ai(£, if) — «|(0, 77)]f/r? + 0(€ 2 ) 

«0 Jo 

(B.53) 


Tlie linearized results given in the text can now 
be derived. Expansion or equations (B.4) and 
(B.5) in powers of e yields 


RTl 

7—1 

rp _J_ Pi L £2 

1 R p 0 y R p 0 Pl 


(B.54) 

(B.55) 


Equation (B.45) is sat isfied identically if a velocity 
potential is defined by the relations 


Differentiation of equation (B.55) with respect 
to t and substitution of equations (B.57) and 
(B.5S) yields 

d7j 1 / dVi 2 ^Vi\ /t> rq's 

7 (Bo9) 

where use has been made of the perfect gas relation 
for the isentropic speed of sound 

<fi=yRT 0 (B.60) 

Substitution of the foregoing relations into equa- 
tion (B.4G) leads to 

(BC1) 

It can be seen that (— =— Qi + 0(e), so that the 
\ a / 1 «o 

last equation can be written as 

(B.G2) 


Ot~ 07) Po 


To obtain equations (l)-(5) of the text from these 
equations two further steps arc required. First, 
the primed quantities in the text are perturbation 
quantities so that, to lowest order, each is € 
times the corresponding quantity with subscript 
1. Since the e would appear as a factor in every 
term, it would cancel. Secondly, if x w (t) is 
neglected in equation (B.53), wo see that x is 

equal to — to lowest order. With this substitu- 
te* 

lion in equations (B.5G)-(B.59) and (B.62), we 
obtain equations (l)-(5) of the text. Alter- 
natively, instead of neglecting x w (t ), we can 
measure x from the wall. That is, we can interpret 
the variable appearing in the text as a new vari- 
able x , with x=x— x w (t), so that £=-- to lowest 

a 0 

order. It will be seen in the following discussion 
that the second alternative is essential in the 




Pi = ~Po 

Substitution of (B.56) into (B.44) .leads to 


<Vi 

(B.56) 

present problem. 

an expansion about €=0 

dy 

Let us consider 

<Vi 
’ dt 


with x as variable 

in place of 77 . As in any 

(B.57) 

acoustic theory, such 

a development would include 


a transfer of the 

boundary conditions from 


dpi 

dt 


— = — p^o 


2 dVi 

C>7? 2 


(B.58) 


x=x w to sr=0; that is, 

du(t,x) 


u(t, x w ) = u{t, 0) + 


dx 


| x *»0 


x„+0(xl) (B.63) 
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It can be seen by integration of equation (B.S) 
that x w is formally of order t. In classical acous- 
tic theory this fact is used to arrive at the con- 
clusion that the boundary condition for u x is 

Wlien x w (t) is large, as in the impulsive-piston 
problem at large t , it is customary to replace x 
with x— x w (f) in the solution. It can be seen from 
equation (B.53) that, to lowest order, this is 
equivalent to using rj as a variable instead of x . 
For the case of a radiating gas, there is an impor- 
tant additional consideration. In appendices C 
and D it is shown that the exact solution (for 
sinusoidal boundary conditions) of the foregoing 
linearized equations is not analytic at the wall. 
It follows that the boundary conditions cannot 
be translated according to equation (B.63), even 
for small x w . Therefore it is essential to use r\ 
as a variable rather than x in this problem. Al- 
though x cannot be used, the variable x — x—x„(f) 
could be, since it would not be necessary to transfer 
the boundary conditions in that case. This is the 
basis for the statement made in the text, that the 
equations obtained by neglecting x w are correct to 
lowest order if x is measured from the wall, rather 
than from a fixed origin. 

One other type of nonuniformity appears in the 


present linearization, as well as in classical acoustic 
theory. A discrepancy appears at large distances 
from the source of a disturbance owing to cumula- 
tive nonlinear effects. This can be illustrated by 
considering the impulsive-piston problem, wherein 
a compression wave travels outward from the wall. 
In the foregoing discussion, it is concluded tlmt, 
in the first approximation, x should be measured 
from the wall. For disturbances in the rear of the 
compression wave this is correct, since they are 
indeed moving into a gas that is at rest relative to 
the wall. On the other hand, the foremost part 
of the wave front is moving into a fixed gas. Thus 
the wave velocity is in error by a small amount at 
the front of the disturbance, an effect which can 
lead to a large error in the relative positions of 
elements of the disturbance in the course of a 
movement over a large distance. It can be seen 
that similar errors will result from evaluating the 
wave speeds at the temperature of the undisturbed 
gas, rather than at the correct local temperature. 
These effects have been treated by a coordinate 
stretching procedure (refs. 25 and 2G). Tn this 
technique, the linearized results are not discarded 
in favor of a new approach. The independent 
variables arc replaced in the solution by new 
variables which arc functions of the original vari- 
ables. It is assumed that this procedure can be 
carried out in the present problem, but the matter 
will not be further investigated here. 



APPENDIX C 

INVESTIGATION OF EXPONENTIAL APPROXIMATION OF ATTENUATION FACTOR 


In reference 1 and in the present work, the 
attenuation factor F{r\) appearing as the kernel of 
the integro-differential equation (30), is approxi- 
mated by an exponential according to equation 
(34). In this appendix the validity of the ap- 
proximation is investigated in the case of a grey 
gas, for which 

(C.l) 

Jo 

It is found that tlie approximation does not yield 
a uniformly valid approximation for the gradient 
of temperature at the wall. The same is true for 
higher derivatives of the other physical quantities. 
However, no nonuniformities are found in the 
approximations for the physical quantities them- 
selves near the wall. At large distance from the 
wall two kinds of error can occur as follows: (1) 
In any acoustic theory, cumulative nonlinear ef- 
fects appear in the evaluation of a flow field far 
from the source of a disturbance, even in the 
lowest approximation. If the disturbance is suf- 
ficiently small at large distances, these effects may 
not be important. When they are important, 
the linearized results must be corrected to obtain 
a uniformly valid approximation (see refs. 25 and 
26). (2) If an approximate solution is used, 

rather than an exact linear result, there may be 
additional cumulative effects which will cause 
error in tlie prediction of the flow field far from the 
source of disturbance. Again, such effects may 
not be important if tlie disturbance attenuates 
sufficiently. These matters will not be investi- 
gated here. In the following, an attempt will be 
made to find any other possible sources of error 
which may result from the exponential approxi- 
mation of tlie attenuation factor. For this pur- 
pose we will concentrate on the solution for the 
response of a grey gas to sinusoidal houndary 
conditions. In that case it is expedient to express 


the potential and tlie boundary conditions in the 
following forms 


PT 

<p- ‘ °Ro [77(y)f> f “' ] 

(C.2) 

CO 


dtp /t>7 'o R< . 

dr Xj=0 Qq 

(0.3) 

Re 

(0.4) 


Re [//(ao)<? f «*] = fmite quantity 

CO 

(C.5) 

where 

y=oix/a 0 (C.G) 

(see ref. 1). The quantities A and B are dimen- 
sionless complex constants, assumed to be speci- 
fied. Tlie magnitudes of A and B must be small 
compared to one, in conformity with the require- 
ment of small disturbances imposed in the 
linearization. 

Substitution of equation (C\2) into equation 
(30) yields 

nty)+H"(y) = -iKb |[/?- 77 ( 0 ) 

-t 77" (0)1 KM- £' EUKy-y)] 

+ '/r' '(?/)] '/.'/+ r E t [b{7,~y) 1 [77' (7/ ) 

77' "(?/)] tfy } (C.7) 

where 

ft' = 8(y— l)<r7 , 2/7?p 0 ao (C.S) 

and 

6=« 0 ^o/ w (C.9) 

The boundary conditions on the complex dimen- 
sionless quantity 771?/) are found from equations 
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(C.2) to (C.5) and arc 

n '{ o)-A (C.io) 

//(»)= finite quantity (0. 11) 

In reference 1 , approximate solutions of equation 

(C.7) are found by replacing the K 2 function with 
an exponential. Here we will consider an approx- 
imation of the form 

E»(y) (C.12) 

i= i 

Solutions can then he found in the form 

(C.13) 

; = 1 

Substitution of equations (C.12) and (C.13) into 
equations (C.7) and (C.IO) and equating the co- 
efficient of each resulting exponential term to zero 
yields the relations 


1+C?- 

*“^0 + 7 )s m, cj-nfi s ° 



for all j 

(C.14) 

J / 

s( 

J=1 v 

fOT “ lu 

(C.15) 



(C.16) 


j=i 


Specializing to the case L— 2, equation (C.14) 
can he written as 


IT Cj— iK0 | 

v ‘ y) tcJ-W 

(i+-)^ 

y 

where 


(C.17) 


7C=2 ~ k 
n i 

(C.18) 

and 



(5~riib 

(C.19) 


Equation (C.17) is a sixth-degree algebraic equa- 
tion in Cj. The six roots arc then the quantities 
Cj in the solution represented by equation (C.13). 
The intcgro-differentiul equation for Tl(y) can be 
converted to a purely differential equation of 
sixtli order. It follows that all solutions of the 


integro-differential equation are contained in the 
general solution, equation (C.13), with J= 6. 
Since the characteristic equation (C.17) contains 
Cj only in the combination c) y half of the roots will 
have positive real parts. The corresponding 
values of C s must be set equal to zero to satisfy 
the boundary condition, equation (C .11). Taking 
the first three roots to be those with negative real 
part, the J in equation (C.13) can then be set- 
equal to three, and equations (C.15) arid (C.16) 
can he solved for the amplitudes C\ y C 2f C 3 in 
terms of the known quantities A, B, c u c 2 , c 3 . 

If mjmi is small compared with one, the roots 
of equation (C.17) can be found by an expansion 
about m 2 /wi equal to zero. The lowest order part 
of the result is that obtained by setting the right 
side of equation (C.17) equal to zero. This is the 
characteristic equation used in reference 1 and in 
the text of the present work. It has previously 
been shown that this equation can he solved to a 
good approximation by an expansion for small 
values of yy— 1. This expansion will be used 
here, since it greatly facilitates the manipulations. 
The results from the double expansion arc 


Ci — —i 


1 + (\'T 1) 


ia 

1 — ia 


m 2 ft- n 

( v7 -d 




[1 — ia] 




/ n 2 fi\ C m 2 iKPlOh/ri,) 2 —]] 

V«i ) \ 2m i («a/«i) 2 — 1 — iKPtruJnx) 2 

+ higher order in (v 7 — 1) and mjm , I (C.22) 
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where 


Kffly 


(C.23) 


7+1 


V 2 


Specializing to the case of constant wall temper- 
ature (jB= 0), equations (C.15) and (C.16) can be 
solved with the aid of the double expansion to 
obtain 

(?!= ta+0( V 7-l)+0 [(V7— 1) (C.24) 

a=0(x 7- 1 ) +0 [( V7- 1 ) ~] (C.25) 
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order, derivatives of C s exp (c 3 ?/) can become large 
at y — 0 when n 2 /S/?h is large. Substitution of 
equations (C.20) to (C.26) into equation (C.13) 
and differentiation three times with respect to y 
leads to the expression 


mi 


'40 

::)■-] 


o -um (' 




+o(^i)+o[(^-i) ( +) 2 ] (CV27) 

The neglected terms are small compared to those 
retained for all values of the parameters. The 
term proportional to ?n 2 /?n 1 is small compared to 
the first for small m 1 lm x ai all values of the 
parameters, except when ?? a /3/?u is large compared 
to one. For that case equation (C. 27) can be 
simplified to 


+ higher order in (v 7 — 1) and m 2 /m l (C.26) 

From these results it can be seen that if m 2 /w, 
is small compared to one, the values of c h c 2f C n 
and C 2 will differ only by small amounts from their 
values at m 2 /m ! = 0 for all values of K and (3. 
Such a change in the wave speeds c-ould cause a 
nonuniformity at large x because it is a cumulative 
effect. For the response of a grey gas to an im- 
pulsive motion of the wall, however, it can be 
shown that such a nonuniformity does not occur. 


Tliis is so because the 


o f-‘( 

Lwi 


Y7— 1) 


corrections 
arge values 


to the wave speeds go to zero at the 
of /? involved in the evaluation of the solution at 
large x. In other words, the component waves 
with appreciable wave speed discrepancies are 
sufficiently damped that their amplitudes become 
negligible at large where their positions arc 
given inaccurately by the single-exponential ap- 
proximation of the attenuation factor. Possible 
nonuniformities of this nature should be considered 
in applications involving nongrey gases or wall 
boundary conditions different than those used here. 

The amplitude <7 3 goes to zero when m a /m i is 
zero and will cause only a small change in II (y) 
when m-tlm-i differs from zero by a small amount. 
However, since c 3 is equal to ( — Wa^/n-j) to lowest 


H' " (0) = —iA— (\7— 1) 5 s A 

m i 
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X-V+... (C.2S) 

Ui 


Using equations (C.2), (2), (3), and (4), it can 
be seen that the first derivative of temperature, 
second derivative of velocity, and third derivative 
of pressure (evaluated at the wall) each have terms 
proportional to 77' "(0). Hence these and higher 
derivatives of T, ?/, and p become singular at the 
wall as (w 2 j3/ni) approaches infinity. However, 
the region of nonuniformity becomes exponentially 
small, since the C 3 term involved is multiplied by 
exp [-(n 2 j3/-n.i)2/]. An effect of this type is pre- 
dicted in reference 1. The basis for the prediction 
can be seen from equation (C.7), in which the 
first term on the right is proportional to E 2 (by). 
The first derivative of the E 2 function is logarith- 
mically singular at a zero value of its argument, 
whereas the first derivative of the integral terms 
of equation (C.7) is not. If follows that the 
only remaining part of the equation, the quantity 
II{y)+II"(y)i must match the singular behavior; 
that is, 77'" (y) must be infinite at y= 0. Further 
details on the form of 1 7"'(y) near y = 0 are 
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given in appendix D, where a procedure for 
obtaining an exact solution is discussed. 

To determine at which values of j3 and y the 
foregoing nonuniformity becomes important, spe- 
cific values of m and n arc needed. Values of 
m u m 2t n Xy and n 2 have been found which satisfy 
the following requirements: 

( 1 ) 

=q=E 2 (Q) = 1.0 

(2) . 

j {m l e~ n ^^m 2 e~ n ^)e~ n ^dri= I E 2 (rj)e~ n ' 2V d7j 
Jo Jo 

(3) 

(mie~ n ^^m 2 e~ 7l ^)rjdr}— I E 2 (y)vdv~o 

Jo Jo o 

( 4 ) 

Least squares fit of m x e to E 2 (rf) 

The first two of these conditions were chosen to 
obtain a good fit of the approximating function 
to E 2 (rf) at small values oT 77. The third require- 
ment insures a correct result in the Rosseland 


limit of strong absorption. The last condition 
is self-explanatory. The resulting values of m 
and n are 

Wi^0.745 
n i = 1.532 
m 2 = 0.255 

71 2 = 20. 

The second term in equation (C.28) has its largest 
value when Kj3 and a are of order 1. For a value 
of y—7/5, and the foregoing values for m and n 9 
the second term begins to exceed the first at about 
/?— 1.0. The region of nonuniformity at that 
point is confined to values of y less than about 
0.1 and is confined to y less than 0.1 //3 for larger 0. 

Three types of possible nonuniformities arising 
from the present approximation procedure have 
been disclosed in the investigations contained in 
this and the preceding appendix. The most 
serious are probably those that can occur at large 
distance from the wall. In the classical acoustic 
theory, nonuniformitics of this type have been 
treated by a coordinate stretching process. The 
details of the procedure for a similar program in 
the present problem are far from obvious. 



APPENDIX D 


EXACT SOLUTION FOR A GREY GAS 


Tn this appendix the validity of the exponential 
approximation of the attenuation factor will he 
further investigated for the case of a grey gas and 
oscillatory boundary conditions. Equations (C.l) 
through (0.7) of appendix O are the exact relations 
for this case. Equation (0.7) is a linear integro- 
diffcrential expiation similar to that appearing in 
the Milne problem (isotropic scattering of radia- 
tion or of slow neutrons; see, e.g., ref. 30). A 
method previously employed for solution of the 
Milne problem could therefore be used here. This 
method, which is exact, utilizes the Fourier trans- 
form plus the Wiencr-TTopf technique for factoring 
the transform. Tn reference 32 a solution based on 
the Wicner-TTopf technique is given for a problem 
even more nearly analogous to the present one 
than the Milne problem. An investigation of 
equation (C.7) by this method would be desirable, 
but- will not be made here. Instead, an alternative 
procedure will be used which yields information on 
the properties of the solution at small and large 
values of y. No numerical examples will be given. 

Tf equation (C.l 3) is substituted into (C.7), and 
the integral expression equation (C.l) substituted 
for the function, it is found that the resulting 
relation can be partially evaluated by an inter- 
change of the order of integration. A solution 
would then be obtained by setting the coefficients 
of each exponential term equal to zero, but there 
is an additional part containing an integral of a 
function of jj . times exp ( — byjpi ) which cannot be 
zero. This suggests tlmt the solution is of the 
form 
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(D.2) 


In evaluating the integrals over jx on the left of 
this equation, the principal part was taken of a 
singular integral. The same must be done in 
evaluating the singular integral on the right side, 
since the two singularities cancel each other. 

If the dummy variable of integration, 0, on the 
left side is changed to ju, equation (D.2) can be 
written as 


^{i +<5 -w(i + f)[ 2 +il.,(^)]} 

[‘2-m In (}^)]} G(n) + -Ub [-7? 


C i e clv 


+g( ,+ ?) r '(i .&>£. 

0 +£)«*> (S>]) 


(D.3) 


n(y)'-=it Cj€ c i v + f 1 G(d)e- bv/9 <l0 (D.l) 

j=l Jo 

When the procedure just described is applied using 
this expression, the result is 


This equation will be satisfied if the following 
relations hold 

+2) [ 2+ i,„(^)]„„ for .11 j 

(D.4) 
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G(p) 


(D.5) 


Equation (D.4) is a transcendental equation 
and lias an infinite number of roots. However, by 
appealing to a principle requiring that the solution 
be a continuous function of the parameters as in 
reference 32, and excluding those roots with posi- 
tive real parts, it is found that only two roots 
survive. The properties of these can be conven- 
iently studied using an expansion for small values 
of y— 1. The results are qualitatively similar 
to those given in reference 1. One of the roots 
leads to a wave speed which differs only slightly 
from that of a classical acoustic wave as in the 
approximate solution. For values of greater 
than one, the properties of the other root are also 
given correctly by the approximation. However, 
the exact characteristic equation (1X4) indicates 
that at values of nb somewhat less than 1.0, the 
second root disappears. The approximate solu- 
tion does not reproduce this effect, but instead 
indicates that the wave speed approaches infinity 
at this point. The approximate result is qualita- 
tively correct in spite of this difference, since the 
integral term in the exact equation (D.l) corre- 
sponds to waves with an infinite wave speed, and 
it is this term that is simulated by the second root 
in the approximate solution for small values of the 
product kI). 

Equation (D.5) is an integral equation for the 
amplitude (7(0) appearing in the integral term of 
equation (D.l). No solution has been found. 
However, if it is assumed that one exists, some of 
the properties of the solution of equation (C.7) 
can be deduced from the form of equation (D.l). 

If an asymptotic* expansion of equation (D.5) 
for small 6 is considered, it can be seen that the 
expansion is of the form 

G(6)~ aJP+brfP In 0- f a^ 3 + . . . (D.G) 

Substitution of this into equation (D.l) shows that 
the third derivative of //(?/) will be singular at 
y = 0 as anticipated in the approximate considera- 
tions of appendix C. 


Substitution of equation (D.l) into the bound- 
ary condition (C.10) yields 

't,cfl,-b\ i yp-d9=A (D.7) 

J=i Jo y 

This relation together with equation (D.5) must 
determine the values of the amplitudes C x and 
Co if a unique solution exists, xSince the function 
G(6) has not been found, exact values of the 
amplitudes C L and C 2 are not available. Although 
the values of these amplitudes arc not known, 
the functional dependence of the solution for 
large y can be found from equation (D.l) in cases 
where the real part of c { or c 2 is greater than — b. 
In those cases the unknown integral term will be 
exponentially small compared to the leading term, 
winch is of the complex exponential typo. The 
wave speed parameter for this functional depend- 
ence at large ?/ is known exactly, since it is a root 
of the characteristic equation (D.4). 

The characteristic equation (D.4) is contained 
in a result given in reference 7. There the equa- 
tion is more complicated, since it includes the 
effects of viscosity, thermal heat conduction, a 
finite velocity of light, and a frequency dependent 
absorption coefficient. No attempt is made in 
reference 7 to describe the source of 1 he disturbance 
or include the effect of a wall. Instead, the 
development is based on conditions far from any 
obstacle. Also the investigation of t lie roots of 
the characteristic equation is confined to a study of 
the one root corresponding to a modified-classical 
wave. 

Since an exact solution of the present problem 
would be obtained if equation (D.5) could be 
solved, it is of interest- to find whether an equation 
of this type has been discussed in the literature. 
As it stands, equation (D.5) contains the two 
unknown constants C x and C 2 . If a solution to 
the complete problem exists and is unique, C t 
and C 2 must be determined by equations (D.5) 
and (D.7). Tt can be seen in equation (D.5) that 
G( 1) must be zero because of the singularity in 
tlie log term. Using this fact, a relation between 
the given constants A and B can be found which 
will result in <^=0. A similar relation leading 
to =0 can also be found. When either of these 
conditions is satisfied, a fundamental solution 
is obtained containing only a (q term or only a 
C 2 term plus the integral term involving (7(0). 
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The integral equation for G(6) in either case no 
longer contains unknown constants. Tor arbi- 
trary values of - 1 and /?, the solution will then he 
a superposition of the two fundamental solutions. 
The integral equations for <7(0) corresponding to 
the two fundamental solutions can he expressed 
in terms of the functions l\ (0) and f t (0) as follows 

(n - s) 

=u&{[l+JJ/ 2 (0)0'70]( 

G(6) is given in terms of the solutions of these 
equations, and the specified constants A and B , 
by the relations 

<»•>« 


MO) 


(-«- ) 
\b + C$) 


where 
a= 


^)mKI>H <D1,) 


(D12) 

The amplitudes (7i and <?, are determined by 

(D - 13) 

a -(^)cr)[ B+ r *"H <d - i4) 

For values of #cb such that c 2 does not exist, the 
correct result is obtained by setting c 2 =0. It 
then follows that / 2 (0) — 0 and (72 = 0. 

Solution of the complete problem, now depends 
on solution of the integral equations (D.S) and 
(D.9) for the unknown functions /i(0) and / 2 (0). 
These two independent equations are essentially 
of the same form. It was pointed out to the 
author by Harvard Lomax, of Ames Research 
Center, that a general solution for equations of 
this type is given in reference 33. However, one 
of the conditions used in the derivation is not 
satisfied by equations ( D.S ) and (D.9), because 
of the singularity at 0—1 in the log terms. There- 
fore, it is not known whether the solution is valid 
in the present case. Alternatively, if it is assumed 
that f(6) is bounded in the interval O<0<1, an 
approximate solution could be obtained by means 
of a truncated expansion of/(0) in Legendre poly- 
nomials. The result from this procedure would 
be of interest for comparison with the results from 
the exponential approximation of the attenuation 
factor used in the text, hut the matter will not be 
pursued here. 



APPENDIX E 


E\ALUATION OF INTEGRALS FOR LIMITING CASES IN THE IMPULSIVE-PISTON PROBLEM 


In the text it was shown that, for the impulsive 
problem, all disturbance quantities are zero at 
slightly positive values of r. The amplitudes of 
steps which occur in the velocity, pressure, and 
temperature were- also evaluated. Using similar 
methods, the variations of the flow quantities at a 
point far from the wall can be found. The 
behavior of the solution for very large values of 
the radiation parameter k can also be found in 
closed form. Finally, a result for very large r can 
be derived. 

For the velocity at a point far from the wall, 
the quantity to be evaluated is 

lim - -^^ "lim Tie -- f {A { e c ^AA 2 e C2V ^)e ivT — 

*-■»» i £->w 7T J 0 “ v 

-lim Re “T P 

t > OT 7T LJo J £- 1/4 

+ |* { (AiC'-H A 1 e c rt)e i ’ r< ^\ (E.l) 

As discussed earlier, the integral near v—0 is 
taken on an infinitesimal quarter circle below the 
origin. We wish to show that the contribution to 
the integral from values of v greater than £~ 1/4 goes 
to zero in the limit as £ goes to infinity. This is 
so because the real parts of c |p£ and c^rfgo to minus 
infinity in the limit, but the proof requires a 
knowledge of the behavior of the real parts of 
and c 2 as functions of v for all Jc. 

For the purpose just stated, it can he shown that 

Rer,^<-yj£v£+ e (!) 

Re c,i>t <— i — +L £ 3 / 4 -j- e (£) 

wliere e(£) is a quantity which goes to zero in the 
limit as £ goes to infinity. It follows that 
40 


forf^ocf 
and k^ 0, k finite 


, v 

|eS-f| <exp ^ — Q- ^+1 | 3/ '-e(£)J 

for 

and k^Ojk finite 

Since A x and A 2 are algebraic functions of p, the 
contribution to the integral from integration in the 
range p — £~ I/4 to v— £ will go to zero exponentially 
in the limit as £ goes to infinity provided that k is 
finite and nonzero. 

The contribution to the integral from values of p 
between £ and oo can be shown to be zero in the 
limit by expansion of the integrand for large v . 
Then equation (E.l) becomes 


lim Re 



(Aie Clv t+ A 2 e r * v i)e i,fT — 


(E.2) 


This integral can be evaluated by expansion of 
the integrand for small p. Expansion of equations 
(S2), (S3), (SS) and (89) for small v yields 

Ci — — i — kv - f 0(p 2 ) 


C2 ^ _ (V) V~iv +0(v7) 

A,=i+o(y?) 

-i- o( v >) 


Substitution of these into equation (E.2) gives 
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The integral of the O(-vF) — term is less than 
0(£ -I/8 ) and hence goes to zero in the limit, 
leaving 


lim 


»0,£) 

U 


=lim Re 

£ — > co 


-i p - I/4 

IT Jo 

“■’HSl)**’’] 


P ivir— 


(h 

V 


By reversing the steps leading to this expression, 
the full interval of integration can be restored and 
the result written as 


lim 

£->ro 


2i S-. 


(E.3) 


The second term can be evaluated by contour 
integration and corresponds to a unit step at 
T = £. In this process the path of integration is 
taken to pass below the origin as discussed earlier. 
The first integral in equation (E.3) can be found 
in reference 34. The final result for the velocity 
at large distance from the wall is 



Equation (E.4) indicates a smooth transition 
from a value of zero at r=0 to a value of 1.0 
when r becomes large compared to £. The transi- 
tion is essentially completed within a range of 


values 


of r from 



to 


T=$ + 


a !( 7 — which is a narrow region compared 

V \7~r 1/ 

to the large distance from the wall under 
consideration. 

In the derivation of equation (E.4), the integral 
of the A 2 term was found to be zero. It is some- 
what surprising that an integral involving only 
the A x term yields a nonzero value of the dis- 
turbance at t — £<0. The results from the solu- 
tion of the oscillating piston problem indicate that 
the maximum possible velocity of the modified 


classical waves is the isentropic speed of sound a Q . 
In the present problem the term represents the 
contribution from the modified-elassical wave sys- 
tem, and in terms of the dimensionless coordinates 
r, £ it follows that the maximum velocity of such 
waves is 1 .0. The region r— £<0 is a part of the 
r, £ plane which cannot be reached by waves 
initiated at t — 0, £=0 and traveling at a maxi- 
mum speed of 1.0. The A 2 term, on the other 
hand, represents the contribution from the radia- 
tion-induced wave system. The oscilla ting-piston 
solution indicates that the maximum speed of 
these waves is the velocity of light (taken to be 
infinite). Hence the region r— £<C 0 can bo 
reached by waves associated with the A 2 term, 
but not by those associated with the A x term. 

The foregoing considerations tend to cast 
doubt on the choice of the real axis in the complex 
v plane as the path of integration for inversion 
of the Fourier transform of i/(r, £). Neverthe- 
less, this choice was shown to correspond to the 
boundary conditions and initial conditions for 
the problem under consideration (in the discussion 
following equation (67)). In that discussion it- 
was noted that a branch cut associated with the 
dependence of c x and c 2 on complex values of v 
could be disregarded as far as the integrand is 
concerned. This is true because the quantity 
Ax exp (cxvQ+A* exp (c 2 vf) is continuous across 
the branch cut in question, which is the one in 
the lower half plane. However, A\ exp (civ£) alone 
is not continuous across this cut, but instead 
interchanges roles with the other term A 2 exp 
(c 2 i>£) in passing from one side of the cut to the 
other. As a result, if we wisli to compute the 
contribution from the modified classical wave 
system alone, the path of integration must be 
altered to include an integration around the branch 
cut in order to insure that the disturbance as- 
sociated with this wave system is zero for negative 
time. An integration around the branch cut 
would also be required to evaluate the contri- 
bution from the radiation-induced wave system 
alone. This additional contribution from the 
A 2 term would not be zero, even though the 
integral of the A 2 term along the real axis is 
zero. Instead, the additional contribution from 
the A 2 term is just equal and opposite to the addi- 
tional contribution from integration of the A x 
term around the branch cut. This last fact 
explains why, in evaluating the total disturbance, 
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the integration around the brunch cut can be 
dispensed with. But it also means that part of 
the total disturbance arises from the A 2 term, 
even though the integral of the A 2 term along the 
real axis is zero. Then it is the radiation-induced 
wave system which is responsible for the disturb- 
ance at r — £<C0, a region which cannot be reached 
by the modified classical waves. 

The corresponding expressions for pressure 
and temperature far from the wall are similar to 
equation (E.4); namely, 


P'(rA) ^ 


r 

i 


1+erf h 
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(E.G) 


The term of order (&r/£ 2 ) in equation (E.G) comes 
from the A 2 part of equation (SO). For very 
large r of order £ 2 /&, this term cancels the other 
part of equation (E.G), and the perturbation 
temperature returns to zero. The velocity and 
pressure do not change in this process. Because 
of the homogeneity of the gas and the presence 
of the radiative heat transfer process, the gas 
must reacli a uniform state at very large time. 
In this state the velocity and temperature of the 
gas must be the same as the velocity and tempera- 
ture of the wall. Since the wall temperature 
is held fixed while the wall is moved impulsively, 
the final perturbation temperature is zero. It 
is interesting to note that the final perturbation 
pressure is the same as that which would occur 
in the absence of radiation, even though the final 
temperature is not. 

Another evaluation can be made in closed form 
as a check on the numerical calculations which 
will be discussed later. At very large values 
of k (high gas temperature) the radiative heat 
transfer will tend to hold the temperature con- 
stant. Under this condition, one would expect 
the disturbance to propagate at the isothermal 
speed ao/VT* In terms of dimensionless coordi- 


nates, the velocity transition would then he 
expected to occur near the line r= yy£. It can 
be shown that this is so for intermediate values 
of £, but not at very small or very large £. For 
this purpose equation (78) can be written as 


lim — p^^lim Re 


k — > 


/v— > 3 




=/»+/*+/» 


(E.7) 


The first integral, which includes integration along 
an infinitesimal quarter circle about the origin, 
can be evaluated in closed form for large k. The 
second integral can also be evaluated l> 3 r means of 
an expansion valid for large k. The remaining 
integral can be shown to he zero in the limit as k 
goes to infinity with £ fixed. 

In disposing of the first and last integrals, a 
knowledge of the properties of the quantities 
Ci and c 2 as functions of r, k } and y is needed. 
Except for an infinitesimal deviation below the 
origin, only real values of v need be considered. 
Values of k from zero to infinity, and values of y 
between 1.0 and 2.0 are of interest. From equa- 
tion (81), the following properties can be deduced 
for the foregoing ranges of v } k , 7 : 

1. The complex quantity Ci has no singularities 
(when the value at is properly defined). 

2. The real part of c { has zeros at and at 
v = 00 . 

3. The imaginary part of Ci has no zeros. 

4. The quantity c 2 has a singularity at v=0. 

5. Both real and imaginary parts of c 2 have 
zeros at v— co. 

G. The real and imaginary parts of Cj and c 2 
are negative or zero. 

By factoring out the singularities and expanding 
about the singularities and zeros, it can be estab- 
lished that 


and 


; [i-«,(,-i) i-ty , 

<_i 

b* 


■yjv (v — ik) 




V v(v—ik ) 


(E.S) 

(E.9) 


where a u b h a 2 , b 2 are positive finite nonzero num- 
bers which can be chosen to be independent of 
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v } k, and 7. To establish that a and b can be 
independent of y and k for values of k approaching 
zero or infinity, use can be made of a power series 
expansion of equation (81) in powers of (7 1). 

The resulting power series converges for the values 
of v, k, and 7 under consideration. 

The content of equations (E.8) and (E.9) can 
be expressed by the relations 


Ci^—i 


[i-kMa+jy (E10) 

, (Ell) 

y — ik) 


where b x and b 2 are independent of v , k , and y. 
Equations (E.10) and (E.ll) can then be used as 
upper or lower bounds with appropriate values of 
b t and b 2 . These expressions for c v and c 2 are of 
the same form as those resulting from a truncated 
expansion for small (7 — 1) used in reference 1. 

The same procedure can be applied to the 
expressions for A\ and A* using equations (E.10) 
and (E.ll) to obtain the relations 


i 1 l ( n ( v V**) /XT' 101 

AjrsI— 63(7— 1) (E ' U) 
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These expressions provide upper bounds for the 
quantity l—A l = A 2 with 63 independent of v } k } 
and 7. 

To evaluate the first integral in equation (E.7), 
equations (E.10)-(E.13) can be expanded for small 
v . The results are 
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Substitution of these into equation (E.7), replacing 
kv with r and expanding for small values of r/k 
and yields 
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Evaluation of this integral, including the in- 
finitesimal quarter circle about the origin, leads to 
the result 


/i = 


I-hliin 0 

" /C-> co 



(E.17) 


Evaluation of the third integral in equation 
(E.7) can be accomplished by expansion of equa- 
tions (E.10)-(E.l3) for large v as follows: 

c,«-t-i 1 (7-l)^t[l+o(J)] (E.1S) 
c 2 ~ — bof^v(v— ik) (E.19) 

.4^1-d^oA) (E.20) 

7 3 =lini Re — - f 4 exp 

/;-) oj TV J * t. 

[-Wr-l);ty+o(i)}? (E.21) 

The integral of the 0 ^ term goes to zero in 

the limit. The remainder can be evaluated by sep- 
arating the argument of the second exponential 
factor into real and imaginary parts. After the 
substitution s=vjk is made, the result is 
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Since &i is of order one, this becomes 


4-lim 0 (E.22) 

For the remaining part of equation (E.7), the 
exact expressions for c and A given in equations 
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(82), (83), (88) and (89) must he used. How- 
ever, an expansion for large k with v fixed can be 
made as follows: 


A +1 +„ 2 \ 

-1 27 1 

h \ h / 


+0 (k’fc) 

(E.23) 

c 2 =o( A3) 
\V kvj 

(E.2 14) 

=l-4 s -o[(l) W ,^] 

(E.25) 

f * [ e c '* 

Jh-i 


V 

(E.26) 


By breaking the integral into two parts, Ijk to 1.0 
plus 1.0 to k, and expanding the exponentials in 
power series, the term containing A 2 is found to be 
at most of order £/&. Similarly, the higher order 
term in equation (E.23) can be shown to lead to 
terms which are at most of order y 1/(7 — 1) 3 &£. 
Then equation (E.26) becomes 

7 2 — lim S Ke — - f 

*->« * J k -l 

[“Ct) £ 3 7+° [!' 

(E.27) 


By reversing the procedure in which the h and / 3 
integrals were removed, the full interval of inte- 
gration can be restored and I 2 written as 


(2 7 J_. e ™ “K [-(13) s l] 7 

+°{|U(7- 1 )»i«]-^}) (E.28) 


This integral is of the same fomTas that evaluated 
in equation (E.3). 


Finally, the velocity at large k is found to be 


lirn 


«(*•»$) 

u 


=lim 

A-><» 


“1 

9 


+0 


^1+erf 

{I- 



The corresponding results for pressure and tem- 
perature are 


lim 


pjizAl 

ptflffj 


lim 

00 


■ i 

2vT 


^1+erf 


{ 2 [(S)!] m }) + ° { 


i 

k' k’ 


[(7-imr 


■}’ 


and 


(E.30) 


lim 

£-><» 


aJJ 


lim 



[(7-l) 3 ^] 


I) 

(E.31) 


Equations (E.29)-(E.31) represent the isothermal 
disturbance propagating at the isothermal signal 
velocity as was anticipated. The error function 
transition occurs in a narrow region with width of 
order £/£, which goes to zero in the limit as k goes 
to infinity with £ fixed at a finite, nonzero value. 
In the derivation of these relations, it was ncces- 
s ary to n eglect terms of order r/k, £/£ and 
\\ f(y-iyk£. In fact, the results given here do 
not apply at values of £ of order k~\ where there 
is a boundary layer. Also at large values of r 
and £ of order k } different results are obtained. 

Expressions valid for small £ were given in 
equations (94)-(9f>). Those results are also valid 
for the present case of large k provided that £ is 
small compared to k~\ The velocity, tempera- 
ture and pressure jumps occurring along the line 
r=£, evaluated in equations (97)— (99) , are still 
correctly given in the case of large k. These 
discontinuities are within the boundary layer 
which develops for large k, since they decay ex- 
ponentially in a distance of order k~ l from the 
wall. Also the results for large £ given in equa- 
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tions (E.4)-(E.6) arc valid in the limit as k goes to 
infinity, with | going to infinity faster. 

The values of velocity, pressure, and tempera- 
ture for very large time can be found by taking 
the limit as r goes to infinity. The factor e tvr then 
oscillates rapidly, and the contribution to each 
integral vanishes for all values of v except those 
near singularities in the remaining factors of each 
ntegrand. The only such singularities arc at 
y=0. Expansion of equations (S2), (83), (88), 


and (89) about r=0 and integration of equations 
(78) -(80) yields 

lim 1 %^=1.0 (E.32) 

T — »co U 

lim2^M?-=1.0 (E.33) 

r-»« PoG'ok 

]ini A ) r(r, j ) =Q (E.34) 

7 — > oo a 0 L 



APPENDIX F 


APPROXIMATE CLOSED-FORM EVALUATION OF THE VELOCITY RESPONSE 
TO AN IMPULSIVE MOTION OF A PISTON 


In reference 1, the response to sinusoidal motion 
of the wall was discussed with the aid of a trun- 
cated expansion for small values of (7— 1). The 
same procedure is useful for obtaining further 
information on the impulsive-motion problem. 
Equations (SS) and (89) can be interpreted as an 
expansion for small values of (7— 1) in addition to 
the previously given interpretations of small r, 
large v, and large k . The expansion converges 
for values of 7 between 1.0 and 2.0. If terms of 
order (7— l) 2 are neglected, equations (SS) and (89) 
become 


-?■ £1— (7— 1) 


ikv 


v * — ikv~\~ 1 

Ci=_ V^Sfc) _ 1+(T_1 h 


-0(7— 1) 2 


ih 




+ 0(7-1) 


■J 


Further expansion of the square root factors 
yields 



These expressions for c 2 and c 2 are identical to 
equations (68) and (69) of reference 1, except for 
the use here of the variables k and v in place of 
K and The equivalence can he established 
with the aid of equations (39), (40), and (76) 
herein, and equation (66) of reference 1, 

Si n c e ^ 7-1- ( Vy + 1 ) ( Vy - 1) = 2( Vy - 1 ) + 

0(v^-l) 2 , the expressions for and c 2 can be 
written alternatively as 



The last forms are chosen instead of the previous 
ones because the resulting truncated expansion is 
exact for >°o as well as for k—> 0 , and 

r— >0. Using equations (F.l) and (F.2) in a cor- 
responding expansion of equations (82) and (S3) 
yields 

1 n _ 4 _ 2 (\y—l ) (v— ?’) yT ( \v — ik — \v) 

7 " (p*-ikv+iy 

+0(v^-l) 2 (F.3) 


The factor 7 in the denominator could be deleted, 
but this would increase the value of the error in 
the truncated expansion by a rather large factor 
for a value of 7—1.40. The coefficients of the 
higher order terms in equations (F.l)- (F.3) 
have maximum values of approximately 1.0 and 
arc zero for limiting values of k and v. 

There is another precedent for an approximation 
of the foregoing type in addition to that afforded 
in reference 1 . The isothermal signal speed is a 0 /V Y 
where a 0 is the iscn tropic speed of sound. Thus 
an expansion for small values of (yy — 1) can be 
interpreted as an expansion for a small fractional 
difference between the isothermal and iscn tropic 
speeds. T11 reference 35 a similar approximation 
is introduced for the case of a chemically relaxing 
gas in the absence of radiation. There the approxi- 
mation is based on the smallness of the fractional 
diffe rence between the frozen and equilibrium 
speeds of sound. Comparison of the results of 
reference 35 with the exact results of references 10 
and 31 shows the approximation to be accurate to 
within a few percent for the cases considered. 

Estimates of the error introduced by this approx- 
imation in the present problem can be obtained 
by repeating the closed-form evaluations for limit- 
ing cases previously made. For the jump condi- 
tions at t=£, the following results are obtained for 
comparison with equations (97)-(99) : 


lim 

0 


[ 


«(£+«, £) 

V 


U 


] 




(F.4) 
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« >o L Po fl / PShU 


= e -(W-i)»«+o(v7-l ) 2 (F.5) 
lim r^L+«- {) f) "| 

« xi L R(/ p u r J 

-- ( 7 ^ 1 ) r-^-»*f+0( V 7-l) 2 (F.O) 

The comparison shows Unit equations (F.4)-(F.6) 
are exact at £=0, but the decay factor in the exact 

equations is exp — rather than 

exp — (77 — 1 )&£ as given by the approximation. 
For 7=7/5, the factor (V 7 +I )/ 2 is equal to 1.092. 
The discrepancy does not lead to a nonuniformity 
in the approximation even at large values of &£, 
since the discontinuities die out at large fc£. For 
y— 7/5 the largest error occurs at k£— 5.46 where 
the approximate decay factor is 0.3GS compared 
to the exact value of 0.403. The maximum error 
is then 0.035, which is 3.5 percent of the total 
velocity transition in the disturbance (equal to the 
dimensionless wall velocity which is 1.0). 

Using the approximate equations (F.1)-(F.3), 
the disturbance at large £ is found to be 


lim 


/l+orfP— 3=L 1\ 


+0(v7— l) 2 (F.7) 


Port 


^i)«i/n-c r rr-.^zun\ 

Poflort 2\ 7 L2 v '(v' 7-1)^J/ 


+0(v'7-l) 2 (F.S) 


ht'( T ,i) 17-1 r iw r T ~t— 1 \ 

h: aJJ 2 — V + pi f L 2 v TV7- mil 


+0(y / 7— l) 2 (F.9) 


Comparison with the exact equations (E.4)-(E.O) 
shows that these results from the truncated expan- 
sion are the same except that the factor \(\ 7 — 1)&£ 
in the denominator of the argument of the error 


function is rei 


/ >7 // in tlie 

7 + 1 


exact equations. For 7 


7/5, the factor - A 


equal to 0.910, representing an error of about 10 
percent. However, again the maximum resulting 
error in the perturbation quantities is a considera- 


bly smaller fraction (2 percent) of the total change 
in the transition. Since approximate expressions 
for the wave speeds were used to derive equations 
(F.7)-(F.9), cumulative errors might have been 
expected to cause large discrepancies at £- >0 ° • 
However, this does not occur because the compo- 
nents with the largest errors in wave speed are 
sufficiently damped that their amplitudes become 
negligible before the error in position becomes 
appreciable. 

The approximate results for the limit as k 
goes to infinity with t and £ fixed arc 


lim 


u(j , £) 

U 


=lim 



r — \7£ 

2 vfv^y -!)£/* 



4 0(y7-l) 2 (F.10) 


lim 

k — ) 00 


£+£+)= 1 ; 
U " 



+0( V 7-1) 2 (F.ll) 


lim /|, ' / ’ / ++l = 0(v+-1) 2 (F.12) 

k->™ rt(jt 

Comparison with the exact equations (E.29)— 
(E.31) shows these results to he correct except 
that the factor V(v7 — 1 )£/& in the denominator 
of the error function is replaced by 



in the exact relations. With 7 — 7 /5, the factor 
(\ 7 +l)/ 2 V 7 is equal to 0.923 representing an 
error of 8 percent in position, but only 2 percent 
in value at a given position. 

The results for the limit as r goes to infinity 
compared with the exact equations (E.32) and 
(E.33) show the approximation to be exact in 
this limit. The approximate counterpart of 
equation (E.34) is correct to 0 (^— 1 ). 

The only remaining exact evaluations available 
for comparison are those for the limit as r goes 
to zero given in equations (94) (96). There it is 
seen that the perturbation velocity, temperature, 
and pressure are zero at r - 0 , except at the wall, 
where there are discontinuities. The truncated 
expansion for small values of (\ 7 — 1 ) yields small 
but nonzero values of the perturbation quantities 
at r^ 0 . This slight inconvenience can be re- 
moved without affecting the results for the other 
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limiting cases as follows: In the derivation of 
the result for small r in the exact case, use was 
made of the fact that the disturbance is zero at 
negative r. The factor e ivT was replaced by 
e ivT —e~ tv \ T \ } the second term leading to zero 
contribution because it corresponds to an evalu- 
ation of the original integral at negative values of 
t , where the disturbance is zero. This alteration 
can also be made in the exact equations (78) -(80), 
and the results will be unchanged for all values 
of r and £. In other words, equations (78)— (SO) 
can be replaced with the relations 


u(t, a) 

i 

u 

2tt 


1 


2tt 



(F.14) 


7lT'(r,t) 1 


a 0 L 


2w 




+ 




dv 


T- e -i»\T\) ill (F.15) 


'When the truncated expansion for small (y'y — 1) 
is used in these equations, the resulting perturba- 
tion quantities are zero at r<0, as they should be. 
It can be shown that the evaluations of the other 
limiting cases using the truncated expansion are 
not affected by this change. Hence the previous 
comparisons between approximate and exact 
results remain valid when equations (F.13)-(F.15) 
are used in place of (78) (SO). 

To obtain a qualitative view of the over-all 
flow field, either further approximation or machine 
computation is necessary. Both methods will be 
used. Only the velocity field will be considered 
in this study. For the machine computations 
one might expect that the integrals could be 
evaluated without resorting to the expansion for 
small values of (yy — 1). However, this is not 
feasible, if a machine program valid for all values 
of the parameters and variables is desired. For 
such computations, considerable knowledge of 
the properties of the integrand are needed for a 
proper design of the integration procedure. Also 
the machine computing time required is not negli- 
gible. Since the truncated expansion for small 


values of (\y— 1) yields qualitatively correct 
results, utilizes a simplified integrand, and de- 
creases the required machine computing time by 
about a factor of ten, it will be used here. 

Using equations (F.l) (F.3), equation (F.13) can 
be written as 


R c (r,g) tt,(r,f) 

U V ^ 0 

where 


o(V^-i) 2 


(F.1G) 


U 

(e ivT ~e~ iv ^ T ^)e~ t ^e~ ( dy~ l ' )k ^ 2 ^ 1 ^ vi ~ ik ^ — (F.l 7) 

V 

2(Vy — 1 ) j 
U 7 2tt 

C w (e tvT — e - f 'M) (y 1 —ik/v— 1 ) 0 i) (e e rt—e e rt) 

J-„ (HV-rtiO 2 dv 

(F. 18 ) 

An integral of the form of equation (F.l 7) cannot 
be evaluated directly by machine computation 
because of the singularity at Also the 

infinite interval of integration cannot be treated 
by machine calculation. Therefore it is necessary 
to subtract terms from the integrand which match 
it at r=0 and at large v. Such terms should at 
the same time be simple enough that they can 
be integrated in closed form. If care is taken in 
the selection of the terms to be subtracted, a 
uniformly valid closed-form approximation may 
be found. 

In equation (F.17), the quantity to be matched 
is the last exponential factor, which can be 

written as 

e~ ( Vt-d****/ a -t m = e - av*/! (v-ib) (»+</&)] (F. 1 9) 
where 

X=(Tfy—l)k£ ( 101 ) 

6=(*y2)+V(*/2)*+l (102) 

An approximation of equation (F.l 9) has been 
found which yields correctly all of the results for 
the limiting cases previously discussed. This 
approximation is 

e -x v 2 /\( v - m ( p+i/b )] — € -x/b 2 ^ € -Xv 2 

+ (e ~ x/i > 2 — e ~ x ) e ~ x ^ b2 e ~ iX * /b -\- e ~ x (F. 20) 
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Equation (F.20) was arrived at by considering tlie 
form of the solution for large k. In this limit, the 
quantity b can be replaced by k . However, it 
can be shown that the approximation remains 
valid for small k if b is used instead of k. By 
adding and subtracting the terms on the right of 
equation (F.20) in the integrand of equation 
(F.17), the following results can be obtained 

R ( r i £) ^ a (r j £) j R.yC 7 *? £) | Q^;~ ^2 (F 21) 

where 


UA(T,t) 

U 


f ( e iVT — e~ tvW] ) e ~ * ■*[ (1 — e ~ x/h2 ) e ~ Xy2 
« / “00 

+ ( e -x/i>*_ e -x) e -xM e -tx,/» +e -x] ( k (F.22) 


p-iplrh I 


{ e -Xv*t [(*-i&K P + i/&>] Q e -Xfb*) e -Xp* 

_ 2 U 5 = 1 ) (Vr^ r i )fr-i> , h 

y (1+r 2 — %kv ) 2 / 

(F.23) 

In equation (F.21) the quantity u n {t } £)fU con- 
tains the remaining part of u c {r , £)/Z7 not included 
in u a (tj £)/Z7, and all of u r (r, £)/U. 

Equation (F.22) can be evaluated in closed form 
using reference 34 to obtain 


^(r,{) i r* r i 
U 


Upon substitution of C\ and Co from equations 
(F.l) and (F.2), the result can be written 

«.vO-,i )= 0 (r< 0) 


{ l»2 (J — j— ^ “ ) 

e ~ x (i+ v 2)2+i 

cos i/f— (e -x/!>2 — c- x )e“ Xl ' 2/!|i! cos (s+y)] 

2(Vr^l) j e _, 1F{ cos (x,„£- 0 7 ) — 

7 C a 

cos (X 2 v£— 0 7 )]\^ (r>0) (F.24) 


(v'7— 1 )/~k( 1 +t> 2 ) 

(l + r 2 ) 2 +£-V 


(F.25a) 


Xi=l + Wy-l)W/[( 1 +v 2 ) 2 +*V] (F .25l>) 


g 2 =Jr+ 1 cos (I tan- 1 /[,(1 +tV) 1/4 ] (F.25 c) 
X 2 = Jli 1 sin (I tan- 1 *) /[„(! +*V) 1/4 ] (F.25d) 


r 3 =(l+A-V)' 


r2= \ /_ 2 ~r r > 


09=4 tan -1 f- 


(F.25c) 


(F.25Q 


(F.25g) 


r 4 =[( 7 - 3 cos 0 2 — 1) 2 +Hj sin 2 0 2 ] 1/2 (F.25h) 

0„=tan _1 fr, sin 0 2 /(l — r 3 cos 0 2 )1 (F.2. r )i) 


+^( 1 +^)+! 

' >0 


+- e~ x I 

^2 ' 


JjV 

< 

1! 

*0 

5^. 

(F.2f)j) 

0 5 — — tail 1 (1/r) 

(F.25k) 

— 2 tan -1 [/^/(l + r 2 )] 

(F.251) 

r s — (1 + f 2 ) 2 +£V 

(F.25m) 

d 7 — ^4 + 05 — #6 

(F.25n) 


In equation (F.23) the symmetry of the integrand Equation (F.24) is now in a form suitable for 

can be used to reduce the interval of integration evaluation by machine computation as explained 

and eliminate the imaginary part of the integrand. in appendix G. 



APPENDIX G 

NUMERICAL EVALUATION OF INTEGRALS 


In this appendix, the results of the numerical 
evaluation of equation (F.24) by means of elec- 
tronic machine computation will be discussed. 
Since the singularities have been removed, the 
only remaining difficulty in the numerical pro- 
cedure is that associated with the infinite interval 
of integration. The difficulty cannot be removed 
by a transformation since there arc a large number 
of zeros of the integrand in the interval. In the 
original integral, there were an infinite number of 
zeros of the integrand before subtraction of a 
function which was evaluated in closed form 
(eq. ( 100 )). As a result it is only necessary to 
integrate equation (F.24) over some finite interval, 
beyond which the contribution is negligible. It 
is not feasible to set an arbitrary large interval of 
integration, because the required computing time 
is not negligible. Therefore the approximate 
required limit of integration was found as a func- 
tion of the parameters k } 7 , and the values of r 
and £. To follow the variations of the integrand 
within this interval it was necessary to break the 
interval into one hundred parts. The calculations 
were then checked hy doubling the interval of 
integration and doubling the number of points 
used. Since a large number of calculations were 
made, some effort toward optimizing the program 
was made. Finally a program was devised which 
would lead to valid evaluations of equation (F.24) 
for values of k between 0.001 and 1000 and 7 
between 1.0 and 2.0. A semiautomatic process for 
choosing the appropriate values of t at which 
calculations should be made for given values of 
k, 7 , and £ was included. 

For all values of the parameters and variables at 
which calculations were made, the part of the 
solution evaluated numerically (cq. (F.24)) was 
small compared to the total variation in u(t, Q/U 
(equal to 1 . 0 ) . Therefore the results from equation 
( 100 ) were used in the text to summarize the 
findings from the numerical investigation. The 
results presented there are rcplotted in figures G, 
7, and S on an expanded scale for comparison with 
the more exact numerical computations. In figure 
50 


6 it can he seen that for small k the approximate 
dimensionless velocity u A fU is correct to within 2 
percent of the total transition. These results 
correspond to the data used in figure 4 . The 
evaluations for intermediate A*, used in figure 3 , 
are compared with numerical values in figure 7 . 
Here there is a maximum error in u A /U of about 12 
percent. This occurs at £=LS19 and t=£. The 
amplitude of the step at this point and at all other 
points is given correctly by u A jTJ, but the approxi- 
mation overestimates the magnitude of the pre- 
cursor to the main part of the disturbance. Figure 
S shows the comparison for large k , corresponding 
to figure 5. The discrepancies are similar to those 
for intermediate k. 

The results from the closed-form approximation 
u A /U (eq. ( 100 )) have been shown to be a good 
approximation. Therefore equation ( 100 ) could be 
used as a basis for a qualitative understanding of 
the effect of the neglected cumulative nonlinear 
terms which would be important at large distances 
from the walk For this purpose, values of the 
temperature near t—% and r=y y£ are needed. 
These can be found from equation ( 100 ) using 
equations ( 2 ) and (4). The results for the tem- 
perature, found in this way, will not be correct in 
the entire flow field, but- should be sufficiently 
accurate in the regions where they arc needed for 
the nonlinear correction. 

*=. 00 ! 


*£=.5458 *£=5.458 *£=54.58 



Figure G. — Comparison of approximate and numerical 
evaluations of velocity response (* = 0.001). 
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T 

Fir. the 7. — Comparison of approximate and numerical 
evaluations of velocity response (* = 3.0). 


WAVES IN A RADIATING GAS 

*=1000 

*£ = 5.458 *£ = 30 £-5.458 £/*=5.458 £/*=30 



Ficure 8. — Comparison of approximate and numerical 
evaluations of velocity response (*=1000). 



APPENDIX H 

RESPONSE TO IMPULSIVE TEMPERATURE VARIATION OF A FIXED WALL 
With boundary conditions 

0 f<0 

^ 0=constant <> 0 (H.l) 

u(t, Q)=-u u (t)=0 (II. 2) 


T’„(t) = 


f: 


— i\v(v — ik)e c 2 v l](e i 


e-t‘\r\) ^+0(Vy-l) 

(IT. 8) 


and for an initial uniform state, equations (58) to 
(64) and (73) to (77) can be used to obtain the 
solution 
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e c * v i — — cvA e ivj( — 
\Ci c 2 J v 


(H.3) 


T (H.4) 
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(H.5) 


1 + 


Ca 


( l+ V t+i A ( j +■ Vfn A Ca 

i(H.6) 

The variables and parameters have the same defi- 
nitions here as in the previous problem (see the 
table of symbols in appendix A). 

Using the approximation scheme discussed in 
appendix F, the integrals can be simplified to the 
following forms 


u( 


r fc)— _® JL f* — 

A) -T 0 a » 2 rJ_L^^+rJ 


(eV*— gV*) (e ivT - 


! ) T+°(V7-1) (IT. 7) 


r '(r,t;) _—; p / ry \ 

© ^7r J_ OT \ V v—ik) 
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where 

Ci = — i+(W— 1) 


(«'"—!) ~H>(Vr-l) (H.9) 




^47l+0(V7-l) 2 (II- 10) 


c 2 =- 


\v(v—ilc) 


+ 0 (y i 7 — 1) (H. 11) 


Equation (H.9) can be evaluated approximately, 
in a form that can be shown to be correct for all 
limiting values of the variables and parameters. 
The result is 


*0— e~ l ) 


1-orf hU) +e~'( 

T 

1 

\2 A ^7/J \ 

/ 


T'{t,Z) 
0 ‘ 


(TI.12) 

Equation (H. 12) indicates that no discontinuities 
in temperature develop. It can be verified from 
equation (H.9) and also from equation (H.5) that 
there is no discontinuity in temperature at r=£, 
where one might be expected. Also the following 
qualitative behavior of the gas temperature can 
be deduced from equation (IT. 12). At a point 
near the wall (small £) the last term dominates. 
This function indicates a rise in perturbation 
temperature from zero (at t= 0) to a value equal 
to the wall temperature, in a time when r becomes 
of order k h This variation includes a nonzero 
initial slope, wliich can be verified in the exact 
equation (H.5). At a point far from the wall 
(£^^1*0) > the first term in equation (H.l 2) domi- 
nates. 1 his function also indicates a rise in 
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perturbation temperature from zero at t= 0 to a 
value equal to the wall temperature at a later 
time. But here the initial slope is zero, and the 
variation occurs in a region near T=k£ 1 2 3 4 5 6 . Thus, at 
large distance from the wall, the temperature 
variation assumes the character of a diffusion 
process. 

No closed-form approximation for the velocity 
lias been found. However, the following proper- 
ties of the solution can be deduced from equation 
(H.7). The velocity disturbance is everywhere 
small compared to that associated with the 
response to an impulsive wall motion involving 
comparable temperature changes. The velocity 
goes to zero at all limiting values of the parameters 
and variables. No discontinuities in velocity 
develop at r=f. The last two findings can be 
verified in the exact equation (H.3). At a point 
near the wall (small £), the velocity rises with a 
nonzero initial slope, but eventually returns to 
zero. Exactly at the wall, the velocity is zero 
at all times according to the boundary condition 
for a fixed wall (eq. (H.2)). At a point somewhat 
removed from the wall (intermediate £), the 
velocity disturbance consists of two parts. One 
of these is associated with the temperature 
disturbance and, hence, leads to a rise in velocity 
near r—k £ 2 . The velocity returns to zero at 
values of r large compared to kg. Also, as £ 
approaches infinity, this component of the velocity 


disturbance goes to zero for all values of t. There 
is a second part of the velocity disturbance at 
intermediate distance from the wall. This com- 
ponent has a peak near t=% and hence represents 
a compression wave resulting from the nonuniform 
heating. This wave travels at a speed between 
the isothermal and isentropic signal speeds. As it 
progresses, it builds up initially, but subsequently 
decays to zero at large distances from the wall. 
This part of the disturbance also goes to zero at 
large r for all values of f. There is a small 
variation in gas temperature associated with this 
compression wave, hut it is of order (\ 7 — 1 ) and 
is not included in the approximate expression for 
the temperature given in equation (H.12). For 
k—0 or &— > oo the total velocity disturbance goes 
to zero everywhere. 

Since the velocity disturbance is small and goes 
to zero at £-->°o in the present problem, the 
cumulative nonlinear effects will be negligible for 
larger disturbances than they would be for the 
impulsive-motion case. Also from a mathemati- 
cal point of view, it is interesting to note that, 
for the present impulsive wall-temperature prob- 
lem, the linear approximation is uniformly valid 
in the limit of a vanishingly small disturbance. 
In this limit, no discontinuities of the shock-wave 
type develop, although discontinuities in the 
derivatives of the flow quantities do occur at 
r— 0 and at r = 
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